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Abstract 



Parametric coupling of optical and mechanical degrees of freedom forms the basis 
of many ultra-sensitive measurements of both force and mechanical displacement. 
An optical cavity with a mechanically compliant boundary enhances the optome- 
chanical interaction, which gives rise to qualitatively new behavior which can mod- 
ify the dynamics of the mechanical motion. As early as 1967, in a pioneering work, 
V. Braginsky analyzed theoretically the role of radiation pressure in the interfer- 
ometric measurement process, but it has remained experimentally unexplored for 
many decades. Here, we use whispering-gallery-mode optical microresonators to 
study these radiation pressure phenomena. Optical microresonators simultaneously 
host optical and mechanical modes, which are systematically analyzed and opti- 
mized to feature ultra-low mechanical dissipation, photon storage times exceeding 
the mechanical oscillation period (i.e. the "resolved-sideband regime") and large 
optomechanical coupling. In this manner, it is demonstrated for the first time that 
dynamical backaction can be employed to cool mechanical modes, i.e., to reduce 
their thermally excited random motion. Utilizing this novel technique together with 
cryogenic pre-cooling of the mechanical oscillator, the phonon occupation of mechan- 
ical radial-breathing modes could be reduced to (n) = 63 ± 20 excitation quanta. 
The corresponding displacement fluctuations are monitored interferometrically with 
a sensitivity at the level of 1 • 10~^^ m/v^Hz, which is below the imprecision at the 
standard quantum limit (SQL). This implies that the readout is already in princi- 
ple sufficient to measure the quantum mechanical zero-point position fluctuations 
of the mechanical mode. Moreover, it is shown that optical measurement techniques 
employed here are operating in a near-ideal manner according to the principles of 
quantum measurement, displaying a backaction-imprecision product close to the 
quantum limit. 
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1 Introduction 

The tails of comets always point away from the sun. Kepler, in the 17*^ cen- 
tury, already conjectured from this observation that sunlight exerts a force 
on particles in the comet tail. Some 250 years later. Maxwell's theory of elec- 
tromagnetic radiation put this conjecture on solid theoretical grounds, yet, 
collecting experimental evidence for the fact that light carries momentum has 
eluded even the most skilled experimentalists of that age. As a famous exam- 
ple, Crookes attempted to construct a radiometer in which the transfer of op- 
tical momentum makes a vane spin, now famous as the "light mill" (figure [T]). 
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Fig. 1. Radiometer used by Crookes around 1877. The pivoted vanes, each silvered 
and blackened on one face, are free to revolve around the central axis. Crookes 
observed that the vanes started to turn when illuminated by a strong light source. 
The origin of the effect was not radiation pressure, it was mediated by the dilute 
gas still present in the globe. In refined experiments by Lebedev in Russia and Hull 
in the US, effects of radiation pressure on mechanical objects were observed at the 
beginning of the 20*^ century. Image copyright Science Museum/SSPL, London, 
reproduced with permission. 

After much debate, however, it was understood that the observed rotation is 



mediated by the dilute gas in which the vane is kept (Woodruff (1968)). True 
radiation-pressure effects, in agreement with Maxwell's predictions, were not 
observed until the beginning of the 20*^ century, in more sophisticated exper- 



iments carried out by Lebedew in Russia (Lebedew (1901)) and Nichols and 
Hull in the United States ( [Nichols and Hull] ( |1901[ |1903a|bD ). 



The situation changed dramatically with the advent of lasers as highly coher- 
ent light sources in the 1970s. It was soon suggested that one could utilize 
the resulting forces to manipulate the motion of dielectric particles and atoms 
or ions in a controlled manner. Pioneering work of Ashkin at Bell Laborato- 
ries demonstrated trapping of dielectric particles using radiation pressure, a 



technique still widely adopted today and known as "optical tweezers" (Ashkin 



( 1970| )). Soon thereafter, proposals to utilize the laser field to cool the motion 
of atoms emerged. Hansch and Schawlow (Hansch and Schawlow ( 1975])) and 
also Wineland and Dehmelt (Wineland and Dehmelt (1975); Wineland and 



Itano ( 1979D ) suggested ways to use radiation pressure to dampen the ran- 
dom motion of atoms or ions, now widely known as "laser cooling". Ashkin 
and Chu demonstrated this principle to trap and cool neutral particles using 
optical forces (Ashkin (1978); Chu et al. ( 1985] )). In the decades to follow, the 
implementation of these ideas lead to what can be considered a true revolution 
in atomic physics, pivotal for discoveries such as Bose-Einstein condensation 
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and the development of the most precise frequency standards available today. 



Among the most spectacular advances of this laser control is the demonstration 



of quantum ground state cooling of a harmonically trapped ion (Diedrich et al. 
( |1989D ; [Monroe et al] ( |1995D ; [Hamann et~aL| ( |1998D ), which has enabled the 
generation of such exotic motional states as Fock or Schrodinger cat states 
dLeibfried et al] ( [20031 )) of a single trapped ion. However, it is a prime example 
of the enigmatic quantum-classical "boundary" that such quantum effects and 
state preparation have never been accomplished with mechanical oscillators. 
It is an interesting, yet little known fact that the ramifications of radiation 
pressure on mechanical oscillators is a subject that even predates the work of 
Ashkin and co-workers. 



Braginsky's group at Moscow state university has played a pioneering role 
in this regard. Embedded eventually in a very general theory of quantum 



measurements (Braginsky et al. (1980); Braginsky and Khalili (1992, 1996); 



Braginsky et al. (2003)), Braginsky considered already in the year 1967 the 



role of radiation pressure in an optical interferometer (Braginskii and Manukin 
( 1967[ )). Since the end mirrors of a interferometer are mechanically compliant 
they are subject to a radiation pressure force. An important impetus for such 
theoretical analysis has been the use of interferometers as detectors for gravi- 
tational waves. The exquisite sensitivity required in these experiments caused 
a strong theoretical interest in understanding the limits of linear continu- 
ous displacement measurements based on parametric transducers. Braginksy 
identified two important consequences of the radiation pressure force that 
photons exert on the interferometer end mirrors, and has developed a com- 
prehensive understanding of optomechanical interactions as they occur in the 
fundamental building block of a gravitational wave observatory shown in fig- 
ure [2j Monochromatic light trapped in a high-finesse cavity (a Fabry- Perot 
resonator in this case) exerts radiation pressure on the massive end mirrors, 
coupling their oscillatory motion to the light. 



As early as 1967 Braginsky and coworkers (Braginskii and Manukin (1967)) 
recognized that radiation pressure can change the dynamics of the mechanical 
degree of freedom, eflFectively adding an optically induced viscous damping to 
the mirror motion. It was also soon understood that this phenomenon (termed 
dynamical backaction (Braginskii and Manukin (1967); praginskii et al. ( [1970[ ); 



(Dykman 


(1978 


); 


Braginsky 



and Vyatchanin (2002)) - which is in essence very similar to the principle of 
laser cooling of atoms, discovered in the decade to follow. While dynamical 
backaction is a classical eflFect, Braginsky and co-workers also showed that 
the quantum nature of light within the interferometer gives rise to an unsur- 
mountable sensitivity limit in displacement measurements of the mirror: For 
sufl&ciently strong input power, quantum fiuctuations of the radiation pres- 
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Fig. 2. Artist's view of a generic optomechanical system: A Fabry-Perot resonator, 
consisting of two mirrors trapping near-resonant monochromatic hght. One of the 
mirrors is mechanically compliant, here it is mounted on a spring and therefore 
constitutes a mechanical harmonic oscillator. 

sure force induces random motion in the mirror, that masks the displacement 
to be detected ( [Braginsky and Khahh| ( |1992D ; |Caves| ( |1980[ |198lD ) by e.g. a 
gravitational wave, an effect now referred to as the quantum backaction of the 
measurement. While the Moscow group devoted many decades in the develop- 
ment of parametric transducers for gravity wave detection, it was not possible 
to observe the dynamical backaction cooling or amplification, though pioneer- 
ing attempts (Braginskii et al. ( 1970D ) were made by the Moscow group in 
1969. 



Interest in optomechanics reemerged in the 1990 in the context of quantum 
optics, resulting in a variety of proposals exploring quantum eflFects in these 
systems, ranging from proposal for quantum non-demolition measurements 
of the light intensity or single quadrature measurements of the mechanical 



displacement ( Braginsky et al.| ( |1977[ [1980|); [Braginsky and Khahh| ([1996|); 
Heidmann et al. ( 1997D ), as well as the generation of narrow-band squeezed 
light using the mechanical oscillator as an effective third-order-nonlinearity 



(Fabre et al. (1994); Mancini and Tombesi (1994)). Later work has studied 



in great detail the possibility to generate non-classical states of motion (Bose 



et al. (1997)), including superposition state s ([Marshall et al. ( 2003D), e ntangled 
states of several oscillators ( Mancini et al.| (2002); Zhang et al. ( |2003[ ); Pinard 
et al. (2005); Bhattacharya et al. ( 2008b[ )), or motional states entangled with 
optical degrees of freedom (Vitali et al. ( [2007 ); Bhattacharya et al. ( 2008a[ )). 
As a special quantum state, it has also been suggested to cool the mechanical 
degree of freedom to its quantum ground state using cooling by dynamical 
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backaction (Wilson-Rae et al. (2007); Marquardt et al. ( 2007D ; [Bhattacharya 
and Meystrel (|2007D; [Genes et al.| (|2008D; [Dantan et al.| (|2008|). 



Observing such phenomena experimentally however, had been severely com- 
phcated by several experimental challenges. First, optomechanical coupling 
by radiation pressure is usually weak — the momentum transfer of a single 
reflected photon changes the velocity of a (free) gram-scale mass by some 
10~^^m/s only. At the same time, the oscillator displacements associated 
with quantum effects are typically on the scale of its zero-point fluctuations 
where m and are its mass and the resonance frequency. 



respectively. For gram-scale oscillators, such as a small end mirror in an in- 
terferometer, this motion is almost ten orders of magnitude smaller than the 
fluctuations of a trapped atom or ion — usually at the sub-attometer scale. Ded- 



icated experi ments (Bocko and Qnofrio (1996); Tittonen et al. (1999); Hadjar 
et al. (1999); Caniard et al. (2007a)) have been approaching such sensitivities, 
but another challenge persists in all room-temperature experiments: Thermal 
noise tends to mask quantum signatures (caused e. g. by quantum backaction) 
as long as the thermal energy /cbT largely exceeds the energy scale hQ^ of the 
motional quantum (/cb is the Boltzmann, and h the reduced Planck constant). 



In recent years, the tremendous progress in micro- and nanofabrication tech- 
nologies has provided novel opportunities to engineer optomechanical devices. 
For example, toroidal whispering-gallery mode optical microresonators were 
shown to exhibit optomechanical coupling to mechanical modes of the struc- 
ture (Kippenberg et al. (2005); Rokhsari et al. ( 2005] )). The high quality factors 
of their optical and mechanical modes allowed observing radiation pressure 



ampliflcation and coherent oscillations for the flrst time in 2005 (Kippenberg 



et al. (2005); Rokhsari et al. ( 2005[ )). A wide variety of other optomechanical 



systems has since emerged, ranging from the nano- to the macroscale. With 
the advent of these new optomechanical systems, interest in observing their 
quantum mechanical properties (Schwab and Roukes (20051)) revived, an ef- 
fort pursued predominantly with electronic transducers up to that point. The 
new research fleld, now widely referred to as cavity (quantum) optomechanics 
( [Kippenberg and Vahala| ( [20071 [20081)), uses a remarkably diverse set of ex- 
perimental systems as depicted in table [T} which summarizes the properties 
of optomechanical systems studied in the last ten years, reflecting also the 
enormous parameter range covered. 



Within this work, we have developed optomechanical systems based on silica 
whispering-gallery mode (WGM) microresonators. Their compactly co-located 
high-quality optical and mechanical modes render microresonators a very fa- 
vorable system for the study of radiation-pressure induced optomechanical 
coupling. Indeed, we show here that basically all experimental challenges for 
the observation of quantum eflFects can be successfully addressed with this 
system, providing a route to ground state cooling of a mechanical oscillator or 
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measurement of its zero point motion. As a key enabling step towards these 
studies, it has been possible for the first time to demonstrate radiation-pressure 
cooling of a mesoscopic mechanical oscillator based on optically induced dy- 
namical backaction (Schliesser et al. (2006)). Adapting, in addition, advanced 
quantum optical and cryogenic techniques to this setting, we show experi- 
mentally that these systems are capable of closely approaching fundamental 
quantum limits — both in terms of the quality of the displacement sensitivity, 
reaching an imprecision below that at the SQL — and the occupation of the 
mechanical oscillator. 
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2 Theory of optomechanical interactions 



2. 1 Classical description and elementary phenomenology 



To begin the discussion of cavity optomechanics, it is useful to review some 
basic underlying physical concepts and simple limiting cases. Consider the 
generic optomechanical system depicted in figure |3](a). The impinging field 
Sin(t) drives the cavity mode amplitude a{t). Part of the boundary of this 
mode — in the simplest case, one of the end mirrors of a Fabry-Perot cavity — is 
free to move, and its displacement is described by the one-dimensional variable 
x{t). Irrespective of the spatial structure of both the optical mode and the 
mechanical displacement pattern, we assume that the displacement x{t) shifts 
the resonance frequency of the optical mode in a linear fashion. 



where uJc is the cavity resonance frequency for x = and 



90 



dc4 

dx 



(1) 



(2) 



is the optomechanical coupling constant. For the cases depicted in figure |3| 
g'o = —uJc/L for a Fabry-Perot cavity of length L, and = —uJc/R for a WGM 
resonator of radius R. 




Sout(t) 




' 

x(t) 



Fig. 3. Schematic of two generic geometries of optomechanical systems, (a) Lin- 
ear Fabry-Perot-type cavity with a movable mirror, pumped through a partially 
transparent mirror, (b) Whispering-gallery mode resonator pumped by evanescent 
coupling to a waveguide. In both cases, an incoming field 5in(t) drives the intracav- 
ity field a{t). The cavity resonance frequency depends on the displacement x{t) of 
a cavity boundary from its equilibrium position. 
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2.1.1 A moving cavity boundary 



We will first discuss the effect of a moving cavity boundary on the opti- 
cal mode, and neglect the backaction (radiation pressure) of the light. For 
a monochromatic pump wave 5ine~^^^^ the resulting equation of motion for 



the intracavity field amplitude reads (following the notation of (Haus (1984))) 



a{t) = (^-i {uj, + gox{t)) - I) a{t) + ^s.^e'^^^' (3) 

Here, a{t) and Si^{t) are the amplitudes of the field in the cavity, and the driv- 
ing field, respectively, normalized such that |a(t)p is the intracavity photon 
number and |5in(t)P the photon fiux impinging on the cavity. Furthermore, 
uj\ is the laser field angular frequency, the cavity energy decay rate, and 
rjc = To/(ro + Tex) G [0 ... 1] describes the degree of overcoupling, tq and Tgx 
being the cavity photon lifetimes due to intrinsic losses of the cavity and cou- 
pling to the taper waveguide, respectively. As one of the simplest cases, we 
analyze the response of the driven cavity to a sinusoidal oscillation in the me- 
chanical degree of freedom. For x{t) = xq sin(f]m^), the solution for the mode 
amplitude reads 

^oo -H^l + ^^m - ^c) + /^/2 



after all transients have decayed within a timescale of ^ (Schliesser et al. 



( 2008bD ). Here, the Jn are the Bessel functions of the first kind and /5 = 



goXo/Q^ is the modulation index. For small amplitudes Xq, so that \P\ <C 1, 
the intracavity field can be approximated to 

a{t) ^ ao{t) + ai{t) + 0{I3^) (5) 



. . _ 9oXo T^^in / e-^("^+^-^^ e-^("^-^-)^ \ 

""'^^ 2 -zA + /^/2V -z(A + f]J + /^/2 -i{A-n^) + K/2 r 

V ' V ' 

ant i- Stokes Stokes 

where A = c^i — c^c is the laser detuning. Evidently, the moving boundary acts 
as a modulator, building up a pair of sidebands ai in the cavity, with weights 
proportional to the cavity Lorentzian evaluated at the frequencies uji + 
and (jU\ — for the upper and lower sideband, respectively. These sidebands, 
schematically shown in figure |4[ are also often referred to as anti-Stokes and 
Stokes sidebands (Kippenberg et al. (2005D). 



The presence of these sidebands corresponds to a modulation of the intracavity 
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(a) 



Sin{t) 



n 



Sout(i) 



a{t) 




x{t) 



(b) 



On 



4 



Frequency 



Fig. 4. Response of a driven cavity to an oscillating end mirror (a), which leads 
ti the generation of anti-Stokes and Stokes sidebands at frequencies cji + Vl^ and 
c^i — fim, weighted by the cavity Lorentzian (b). 



photon number and it is easy to show that 



(8) 



A2 + (k/2)2 



(1+ 



+ 



A -a 



(A + Q J2 + (^/2)2 (A-l]J2 + (^/2)2^ 
JA + Q J2 + (^/2)2 - (A-l]J2 + (^/2)2^ 



sin(r^m^)+ 



cos(Qm^) I • 



(9) 



The intracavity optical photon number is modulated at the oscillation fre- 
quency f^m, however, the modulation does not necessarily occur in phase with 
the mechanical oscillation as is evident from the cosine-term. In fact, the phase 
lag depends in a non-trivial manner on the laser detuning A, oscillation fre- 
quency and the cavity buildup time Kj , a simple calculation yields for 
the phase lag 



6iag = arg ((?oA(A2 + («:/2)2 - - i^^K)) 



(10) 



which is shown in figure [5} The quadrature component (oc cos(f]m^)) can 
become significant if the cavity buildup time is comparable or larger than 
the oscillation period leading to a viscous radiation pressure force as 

detailed below. Both in-phase and quadrature components are shown in the 
parametric plot in figure [6) 



2.1.2 Radiation-pressure backaction: static phenomena 

So far, only the effect of the mechanical on the optical degree of freedom has 
been considered. However, the richness of optomechanical phenomena arises 
only when the mutual coupling of optical and mechanical degrees of freedom 
is taken into account. The physical origin of the "back-action" of light on the 
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Detuning A//^ 



Fig. 5. Phase lag (/)iag between the oscihation of the mechanical degree of freedom 
and the stored optical energy (or circulating power) in the cavity, assuming go < 0. 
For ±A = ^ ^ the phase lag is exactly =F7r/2. The phase jump at A = is due 
to a zero crossing of the oscillation amplitude of the stored optical energy. 




Fig. 6. In-phase and quadrature component of the intracavity energy (or circulating 
power) modulation when the detuning is varied (colored curves). Different colors 
represent different choices of the cavity linewidth which was changed from 
(green) to 10 Q^n (red) in logarithmic steps. In the blue curves, a fixed detuning was 
chosen (as indicated by the labels), and the cavity linewidth was varied from 
to lOnrn- 

movable cavity boundary is due to radiation pressure. In the simple case of a 
Fabry-Perot cavity, the force arises from the momentum flips of the photons 
reflected from the movable mirror, and thus 

FrAt) = Ht)\':^2hk = = -hgo\a{t)\', (11) 

where hk with k = ui/c is the momentum of the photons. 



12 



The relation -Frp(^) = — %o|fl^(^)P is of general validity and also applies to the 



WGM resonators, as will be derived in section |3.3[ If the cavity boundary is 
free to move, the coupled equations describing the optomechanical system will 
read (in a frame rotating at the laser frequency) 



a{t) = {i{A - gox{t)) - k/2) a{t) + ^Jri^s;n{t) 

\a{t)\^ 



x{t) + T^x{t) + ^i^x{t) = -%o- 



rries 



(12) 
(13) 



where, for the mechanical osciUation, a resonance frequency f^m, viscous damp- 



ing at a rate Fm, and an effective mass meff (cf. section 3.3) are assumed. For 
a constant drive amphtude 5in, these coupled nonlinear equations can be ana- 
lyzed in a first step by finding stable solutions a{t) = a, x{t) = x in which all 
time derivatives (d(t), vanish, requiring simultaneously 



1 



a = 



-i{A - gox) + 1^/2 



and 



(14) 
(15) 



Equations (14) and (15) can be both understood as functions mapping the 
displacement x to an intracavity photon number |ap, as shown in figure [7| 
The self-consistent, physically possible solutions are given by the intersections 
of the two curves. ^ Evidently, the system has at least one solution for arbi- 
trary parameters. For sufficiently high power or finesse, and/or floppy enough 
mechanical oscillators, two additional solutions are physically possible. The 
condition for their appearance is given by 



> 



V3 l^^meff /^^ 
9 ^c%o 



(16) 



which is derived in a straightforward manner from the requirement that the 
maximum slope 3\/3|5inp7yc|5'o|/^^ of the Lorentzian square modulus of (14) 
must exceed the slope mQffQ'^/h\go\ corresponding to (15). Above this thresh- 
old, the optomechanical system displays a well-known bistable behavior, re- 
sulting, for example, in a hysteretic transmission behavior upon the variation 
of the input power. This effect was observed in a pioneering experiment ([Dorsel 



et al. (1983)) at the Max-Planck-Institute of Quantum Optics (MPQ) as early 



as 1983 and is referred to as optical bist ability. Reports in the microwave 
domain followed soon thereafter (Gozzini et al. ( 1985] )). 



Note that the phase of the complex entity a can always be adjusted to fulfill (14), 



as it does not affect (15). 
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(a) (b) (c) 




Detuning laser - empty cavity (MHz) Detuning laser - empty cavity (MHz) 



Fig. 7. (a) Graphical representation of the self-consistent solutions of the coupled 



equations ( [T4| ) (red) and (15) (blue) for the intracavity photon number |ap and the 
radiation-pressure induced displacement x of the cavity radius. Intersections of the 
curves indicate possible stable (full circles) and unstable (empty circles) solutions, 
(b) Plot of numerical solutions for the power |ap/Trt circulating in the cavity and 
the radiation-pressure induced displacement for typical parameters of a silica mi- 
crotoroidal optomechanical oscillator (R = 25 /xm, k = 27r8MHz, = 27r50MHz, 
meff = 20 ng, (jjc = 27r380THz, l^inp = ImW). The dashed line indicates the un- 
stable solutions, (c) The actual detuning A — gox as a function of the detuning A 
of the laser from the undriven cavity resonance. 



2.1.3 Radiation-pressure backaction: dynamical effects 



Quantitatively novel behavior occurs when analyzing the dynamical response 
of the system around an equilibrium (a,x). As first pointed out by Bra- 



ginsky and co-workers in 1967 (Braginskii and Manukin (1967); Braginsky 



et al. (2001); Braginsky and Vyatchanin (2002)), the dynamics of fluctuations 
around the equilibrium not only display new physical effects, but are also of 
experimental relevance, in particular in the most sensitive gravitational wave 
interferometers. Keeping with the simple illustrative approach of this section, 
let us assume the system is in a stable equilibrium if a = a and x = x and an- 
alyze the dynamics of small excursions Sa{t) and Sx{t) if a small external force 
SF{t) is applied to the mechanical oscillator. Substituting a{t) = d + 6a{t)^ 



x{t) = X + 5x{t) into (12)-(13) and introducing the equilibrium detuning 



A = cji - {uJc + 9ox), 
one finds the linearized equations 



(17) 



Sa{t) = (+«A - K/2)Sa{t) - igoa 5x{t) (18) 
eff {Sx{t) + r„ 5x{t) + n'i Sx{t)) = -hgoa{5a{t) + Sa*{t)) + 5F{t), (19) 



where ( [l4j ) and ( [l5j ) were used, and second-order terms oc Sa{t)Sx{t) or oc 
|(Ja(t)p were dropped, as we assume \5a\ <C \a\. Furthermore, without loss of 
generahty, we have assumed real a = d*, which can always be attained by 
adjusting the (physically irrelevant) phase of the incoming driving wave Sin- 
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This equation system is most easily solved in the frequency domain, by apply- 
ing a Fourier transformation to all involved time-dependent variables. ^ One 
then obtains 

-inSa{n) = {+iA - K/2)Sa{n) - igodSx{n) (20) 
-iQ6a%n) = {-lA- i^/2)6a%n)+igod6x{n) (21) 

meff [-n^ - iT^n + Sx{n) = -hg^a {Sa{n) + Sa%n)) + SF{n), (22) 

where 5a* (f]) = {5a{—Q)y was used. Analogous to the previous section which 
considered a moving boundary, we now find that a non-zero displacement am- 
plitude 6x{Q) at Fourier frequency Q induces anti-Stokes and Stokes sidebands 
of amplitudes 

respectively. As a result, the intracavity energy is modulated, giving rise to 
an oscillating force of magnitude 

5Frp{n) = -hgoa {5a{n) + 5a* (Q)) = (25) 

_ _ 2 2 f ^ + ^ A-0 \ . . 

(26) 

The real and imaginary parts of the radiation pressure force in this represen- 
tation are identified as being due to the in-phase and quadrature modulation 
of the circulating power in the cavity. The additional force acting on the me- 
chanical oscillator changes its dynamical behavior, in particular its response 
to the external perturbation. This effect is known as dynamical backaction 
(Braginskii and Manukin ( 1967D ). 



Specifically, substituting (25) back into (22), a modified response of the oscil- 



lator to an external force is found, 

6x{Q) = Xei^{^)SF{Q) (27) 



^ We choose the convention /(fi) = f{t)e'^^^^dt. 
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with the effective susceptibihty Xeff(^), 



fedba(^) 

rrieti 



(28) 



The damping and spring constant induced by dynamical backaction are given 
by 



Tdba — 
^dba = 



k/2 



k/2 



figW 

meffQ V + + («/2)^ (A - + («/2)2 ) 
hgla' / A + Vl A-Q \ 

\{A + + {K/2f ^ (A - ^If + {K/2f ) ■ 



(29) 
(30) 



If the induced changes of the mechanical oscillator's dynamics are small, the 
oscillator still behaves as a damped harmonic oscillator with effective damping 
and resonance frequency (Kippenberg et al. (2005); Schliesser et al. ( 2006| )) 



2^2 



k/2 



k/2 



^gja _^ 

meffJln, \{A + + {K/2f (A - r2^)2 + (k/2)2^ 

hgla^ ( A + Vt^ A-n^ 



2m,sn^ \{A + + (k/2)2 (A - + {K/2f ^ 



(31) 
(32) 



2.2 Formal framework: quantum Langevin equations 



A more general formulation of optomechanical interactions than the simple 
but illustrative considerations in the previous section is possible within the 
framework of a quantum Langevin approach. This enables the full description 
of the quantum dynamics of an optomechanical system. In particular, effects 
related to the quantum nature of light can be treated in an adequate manner. 



2.2.1 Hamiltonian of cavity optomechanics 



Starting point of the analysis is a Hamiltonian formulation of a generic op- 
tomechanical system (Law ( 1995D ). If the mechanical oscillation frequency is 
much smaller than the free spectral range of the cavity, such that a only a sin- 
gle optical mode has to be considered, the system Hamiltonian can be written 
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as 



H = i/mech + ^opt + ^int + ^drive (33) 

^mech = 7, ^ -rrieff^l^x^ (34) 

Hopt = hujc (^a^a + (35) 

^int = %o£ (36) 

i^drive = th ^ {s,^a)e-"^'' - sl^ae+'^'') (37) 



where x and p are the mechanical displacement and momentum operators, and 
d) and a are the creation and annihilation operators of the considered optical 
mode, i. e. n = d)a is the intracavity photon operator, and correspondingly, 
the drive amplitude is now normalized to photon flux at the input of the 
cavity l^inp = P\n/^\' Evidently, this Hamiltonian reproduces the optical 
resonance frequency shift upon mechanical displacement, as 

^opt + Ant = K^c + gox) d)a, (38) 
and simultaneously describes the radiation pressure force, with 

^rp = -^5^ = -hgoa^d. (39) 



2.2.2 Quantum Langevin equations 



From the Hamiltonian, the time evolution of the operators of interest can 



be derived. In addition to the conservative dynamics described by (33), dis- 



sipation of both the optical and mechanical modes, and the corresponding 
fluctuations are taken into account by the following set of quantum Langevin 
equations (QLEs) ( [Giovannetti and Vitali ( 2001[ )) (in a frame rotating at (jJ\) 



d 
~dt 



d{t) 



, -A • -m-m , ('5in + ^'5in(^)) , 5s^^^{t) 

- - a{t) - igox{t)a{t) H — \ (40) 

zy V^ex v^o 



m 



d 



p{t) = -rrieAHt) - hgod\t)d{t) - T^p{t) + 8F,^{t) 



(41) 
(42) 



where the noise terms 5s\^^ ^s^s.ci and 5Fth were introduced. They fulflU the 
commutation relations (Gardiner and ZoUer ( 2004D ) 



[5hn{t),54^{t')] = [5rs,ut),srsl^{t')] = s{t-t') 



(43) 
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and the quantities 



{SsMSslit')) = {SsUt)SsUt')) = 5it - t') (44) 

are the only non-zero correlation functions for the quantum vacuum entering 
the cavity from its two ports: the one through which it is pumped (^^in), 
and the second port (5svac) representing all other loss channels. Here, zero 
thermal excitation of the optical mode has been assumed which is valid in 
the optical domain. In order to adequately describe the mirror undergoing 
Brownian motion the correlation function of the mechanical driving term can 
be shown (Giovannetti and Vitali ( 2001D ) to have the form 



m^{tm^{t')) = hm^sT^ I e-^(*-*')fi Lth + l] g (45) 



As in the previous subsection, the QLEs are simplified in the first place by con- 
sidering static and the dynamical effects separately. To this end, the unitary 
transformations a{t) = d+6d{t) and x{t) = x+6x{t) with {Sa{t)) = {Sx{t)) = 
yields again the requirements ([T4|) ([T5l) for the steady state values of intra- 
cavity field amplitude d and displacement x. 



2.2.3 Dynamics of the fluctuations 



Assuming (a, x) to be known as a stable solution of the system (which can. 



for exa mple, be tested for using the Routh-Hurwitz criterion, see ( |Fabre et al. 



(1994))), the Heisenberg equation of motion for the fiuctuations 5a, and 



5x can be derived. Choosing again the phase of the input field such that a 
is real and positive, and assuming again a strong coherent drive 

d > 1 (46) 

it is possible to derive linearized quantum Langevin equations for the fiuctu- 
ations by dropping terms (x 8a8x^ 8d)8x or 8d)8a., yielding 

4 5a{t) = (+iA- -\ 5a{t) - igoa5x{t) + %0 + (47) 
at V 2/ Jt^ Jt^ 



4 5aHt) = (-iA - -) Sa\t) + igoaSx{t) + ^-^B + (48) 
at \ 2 J Jt^ Jt^ 



^ Sxit) + T^^ Sxit) + niSxit) = -^a (Sait) + Sa\t)) + (49) 
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where the Hermitian property Sx(t) = was used. This set of equations 

is most easily solved in the Fourier domain: 



(-«(A + Q) + k/2J Sa{n) = -igoaSx{n) + 



and 



(5Fth(f^)5Fth[n']) = 27r(5(n + n')hmesTn^n coth 



+ 1 



(50) 



(51) 



meff (fll, -fl'^ - iTra^) Sx{fl) = -%o« {sd{fl) + 5a) (fl)) + 5Fth{^). (52) 
In the frequency domain, 



(53) 
(54) 



(55) 



are the only non-zero correlation functions. Together with the input-output 
relations for the fluctuations 



(56) 
(57) 



these equations constitute the theoretical description of the most important 
effects in cavity optomechanics. 



3 Whispering gallery-mode microresonators as optomechanical sys- 
tems 



Optical microcavities are dielectric structures that confine light to small vol- 
umes for extended amounts of time. In 1989 Braginsky and co-workers in 
Moscow discovered that small dielectric microspheres from glass exhibit opti- 
cal whispering gallery modes that exhibit giant photon lifetimes, correspond- 
ing to ultra high Q-factors exceeding 10^ (Braginsky et al. (1989); CoUot et al. 



(1993); Vernooy et al. (1998bj)). The ultra high Q of microspheres were sub- 
sequently investigated for a variety of applications, including laser stabiliza- 
tion dVassihev et al.| ( 1998])), ca vity Quantum Electrodynamics ( cQED) ([Ver- 



nooy et al.| ( |1998aD ; |Aoki et al.| ( [20061 )), narrow linewidth lasers rtSandoghdar 



et al. 



^(1996)), nonlinear optics at ultra low light powers ([Chang and Campillo 

([1996[ ); [Spillane et al.[ ( [2002[ )) and biophysical sensing (iVo Umer et al.[ ( [2002[ 
2009[ )). That optical microresonators also exhibit mechanical modes that are 
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Toroidal Microresonators 



Microdisk Resonator 



Microsphere Resonator 



Spoke-Microresonators 




Fig. 8. Optomechanical systems based on WGM optical microresonators. Top row 
shows scanning electron micrographs of four different systems studied experimen- 
tally. The corresponding panels in the bottom row depict the displacement patterns 
of mechanical modes which couple to the optical modes particularly strongly. The 
mechanical modes are simulated using the finite-element method (first three panels: 
2-dimensional cross-section, last panel: full 3-dimensional simulation). 



co-located in the same structure and are coupled via radiation pressure, was 



only realized in 2005, in toroidal microresonators (Kippenberg et al. (2005); 



Rokhsari et al. ( 2005[ )). It became immediately evident that optomechani- 
cal coupling in this manner pertains to virtually any optical microresonator 
(Kippenberg and Vahala ( 2007] )) and has been demonstrated in a variety of 
microresonator geometries, including microspheres (Ma et al. (2007)) and mi- 
crodisks or optimized optomechanical resonators ( Schliesser et al.| ( |2008b ); 
Anetsberger et al. ( 2008[ )) as shown in figure [8| In this section, prior to dis- 
cussing the optomechanical phenomena that can occur in optical whispering 
gallery mode microresonators, we analyze the mechanical modes of WGM 
microresonators, discuss their mechanical eigenfrequencies and address the 
question of the sources of mechanical dissipation and other relevant optome- 
chanical properties. 



3. 1 Optical properties of WGM silica microresonators 



An essential precondition to render the weak effects of radiation pressure ex- 
perimentally accessible is a high optical finesse. In a cavity of finesse J^, the 
circulating power — and thus also the radiation pressure force — is enhanced by 
a factor of ^ J^/tt compared to the launched power. Silica microresonators can 
achieve finesses on the order of 10^, exceeding even the best results achieved 
with Fabry-Perot cavities for cavity QED (Vernooy et al. ( 1998bD ). 
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At the same time, in order to reveal dynamical effects of radiation-pressure 
coupling, the photon storage time should be on the order of — or ideally exceed — 
the period of the mechanical oscillator coupled to the resonator. In spite of 
the very short round-trip time of about 1 ps, the storage time can amount 
to several hundreds of nanoseconds owing to the high finesse. This is about 
10 times longer than the oscillation period of the mechanical modes in these 
structures. To ensure well-controlled optomechanical interaction, the optical 
mode spectrum must be well understood. 



3.2 Optomechanical coupling in silica WGM resonators 



3.2.1 Silica microspheres 

The mechanical degree (s) of freedom coupled parametrically to the WGMs in 
silica microresonators are given simply by the intrinsic acoustic modes of the 
structure. Each of the eigenmodes can be viewed, to a good approximation, 
as a damped harmonic oscillator, driven by thermal forces. In the following, 
the nature of these modes, in particular, their eigenfrequency, the damping 
mechanisms, (effective) mass, mode shapes, and coupling to the optical degrees 
of freedom will be discussed for two canonical systems; silica microspheres and 
silica microtoroids. 



3.2.2 Acoustic modes in silica microresonators 

The deformation induced by mechanical modes is described by a vector field 
u{r^t)^ which denotes the displacement of an (infinitesimally small) cubic vol- 
ume element at position r and time t from its initial position. In an isotropic 
homogenous medium, to which no external forces are applied, the equation of 
motion for the displacement field reads ( [Landau and Lifshitz ( 1970D ) 



pU{f, t) = (A + /i) V(V • u{f, t)) + iiV^u{r, t) (58) 
where the density p and the Lame constants 

with a Poisson's ratio and E Young's modulus, characterize the elastic prop- 
erties of the material. 

While for an infinitely extended medium, a continuum of solutions for equa- 
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tion (58) are obtained, for a finite-size body such as a silica sphere or toroid, 
the boundary conditions lead to a discrete spectrum of solutions, such that 
the total displacement 'u(r , t) can be decomposed into modes oscillating har- 
monically at a set of frequencies 



uir 



(61) 



where c^(t) is the displacement amplitude of a mode with index n, Qn its 
eigenfrequency and t^(r) is the spatial displacement pattern of the mode, 
normalized such that 

/y^(fX,(f)rf^r 



(62) 



For more sophisticated geometries (e. g. a silica toroid supported by a silicon 
pillar), it is difficult to obtain analytical solutions for the mode shapes and 
frequencies. ^ In this case, equation (58) with suitable boundary conditions is 



solved using the finite-element method (FEM). As an illustrative example, and 
to validate the accuracy of the FEM results, an analytical solution available 



( 


1906 


); 


Ma et al. 



spherical coordinates, equation (58) is solved in a homogeneous medium by 
functions of the form 



u{f, t) = V(/)o(f, t) + V X $i(f, t) + V X V X $2(r, t) 



(63) 



derived from a scalar potential 00 and two vector potentials $i = (r(/)i,0,0) 
and $2 = {r(p2, 0, 0) with 



'a 



nlm 



(64) 



where g = 0, 1, 2; ji denotes the spherical Bessel function, Yl^ is the spherical 
harmonic function, Vq = ^(A + 2/x)/p is the longitudinal sound velocity, and 

Vi = V2 = \Jl^/p is the transverse sound velocity. The acoustic modes are 
characterized by an angular momentum mode number / (/ = 0, 1, 2, . . .), an 
azimuthal mode number m {—I < m < I) and a radial mode number n {n = 
1,2,...). Here, n = 1 corresponds to the surface mode, n > 2 to inner modes 
and ^nim denotes the frequency of the vibration characterized by the mode 
numbers (nj^m). 

We focus now on the fundamental spheroidal mode (nj^m) = (1,0,0), for 



higher-order modes, cf. (Ma et al. (2007)) and references therein. In particular. 



For simple cylinders, approximate solutions have been developed (Hutchinson 



(1979, 


1980 


); 


Tamura 


(2009)) 



22 



the displacement vector field is purely radial, 



-0 /-X A sin(/ci,o,or) - /ci,o,orcos(/ci,o,or) ^ 

%,0,o(^) = ^0,1,0,0 —2 (65) 

where /ci^o,o = ^i,o,o/'^0 5 ^o,i,o,o is a normalization constant, and the radial 
unit vector. In the following, we drop the mode index (1, 0, 0) for better legi- 
bility, and an index to a vector field now denotes one of its components in a 
given coordinate system. 

To determine the allowed values of k (and therefore Vt)^ the strain tensor 

is evaluated in a spherical coordinate system. For the considered mode, all 
non-diagonal elements vanish, and the diagonal elements read^ 

dur , . (/c^r^ — 2) sin(/cr) + 2/cr cos(A:r) , , 

= ^ = ^ -z (67) 

e,, = ^ = c{t) A ^Hkr)-krcos{kr) ^ ^^^^ 

where c{t) is the excitation amplitude of the mode. The stress is related to 
the strain by Hooke's law, which in a lossless isotropic medium is given by 

Gij = 2iieij + A5ij ^ Skk (70) 
k 

rendering in this case the stress tensor cFij diagonal as well, 

((A + 2/x)/c^r^ — 4/x) sin(/cr) + 4/x/cr cos(/cr) 
Grr = c[t) A — (71) 

{Xk'^r'^ + 2/x) sin(/cr) — 2/x/cr cos(/cr) 
cree = c{t) A — (72) 

{Xk'^r'^ + 2/x) sin(/cr) — 2/x/cr cos(/cr) 
= c(i) A . (73) 

Applying the boundary conditions for a freely oscillating sphere of radius R 
leads to the requirement that the stress is zero at its boundary (r = i?). 



For the definitions of the strain tensor in spherical coordinates see (Landau and 



Lifshitz (1970)). 
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leading to the characteristic equation 



/ _lz|^,^,\tanp_^^Q (74) 
\ 4:vl J kR 

For the parameters of fused sihca, this equation is solved for kR ^ 2.4005 . . . 
or 

^ 2280 m/s 

2^ i^' ^^^^ 

corresponding to a resonance frequency of 91.2 MHz for a 50- /xm diameter 
sphere. Finite element modeling yields the same result within less than 1%, 
testifying to the integrity of our employed numerical simulation. 

Another interesting physical entity is the potential energy stored in the defor- 
mation. It given by 

The diagonal form of strain and stress tensors for a sphere again facilitates 
the analytical evaluation of this integral simple yielding 

U ^ 8.69 • 10^^ -^'R'X^ (77) 
m*^ 

where x = u{R^9^Lp) • = c{t) A{8'm{kR) — kRcos{kR)) / R^ is the radial 
displacement of the boundary, with A ^ 0A27R^ for silica. Figure [9] shows 
displacement, strain, stress and strain energy density of a silica sphere when 
the fundamental mode is excited. Results from analytic calculations and the 
FEM agree very well. 

Figure [TO] compares experimentally measured resonance frequencies of several 
silica microspheres with diameters between 30 and 100 fim with analytical 
calculations. Again, very good agreement is obtained. 



3.2.3 Silica microtoroids 

Toroidal silica resonators, attached a silicon pillar, have fewer symmetries 
than spheres. This structure must be parameterized by several parameters, at 
least by the major and minor radii of the silica torus, the radius of the silicon 
pillar, the thickness of the silica disk, and the offset of the symmetry planes 
of the disk and torus along the z-axis (|Kippenberg et al. (2005|)). The shape 



of the silicon pillar is assumed to be rotationally symmetric, and to constitute 
a quarter circle in the r-z-plane. The radius of this circular arch is another 
degree of freedom, but is usually assumed to be similar to the difference of 



toroid and pillar radii (figure 11) due to the isotropic nature of the silicon 
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Fig. 9. The {n,l,m) = (1,0,0) mode of a silica sphere. Left panels show (from 
top to bottom) the radial displacement u(r^9^(f) • e^, the strains Srr (green) and 
sqq = £(^(^ (blue), stresses arr (green) and aoo = a<^<^ (blue), and the strain energy 
density ^X^zj^^j^^j- Symbols are results of finite-element modeling and lines are 
derived from the analytical calculations, showing excellent agreement. The given 
magnitudes correspond to a 50 /im-diameter sphere containing a total strain energy 
of A:b(300K)/2 . The right panel shows the exaggerated displacement profile (origi- 
nal sphere outlined in black), and magnitude (color coded) as well as the principal 
stresses (indicated by arrow lengths). 

etch. 



The complex boundary conditions (all surfaces free, except for the bottom 
of the silicon pillar) render an analytical solution prohibitively difficult. We 
therefore use a commercial finite-element software to solve equation (58) on a 
discrete mesh consisting typically of more than 10,000 nodes which are auto- 
matically distributed in the simulation volume. Simulations can be run both 
assuming rotational symmetry for the modes, and in full three dimensions. Fig- 
ure [12] shows the obtained displacement patterns for the 19 lowest-frequency 



modes of a toroid of major radius 23.0 /xm, pillar radius 13.23 /xm, minor radius 
2.63 /im, disk thickness 2 /xm and no offset, in a three-dimensional simulation. 
As shown in figure 13, the simulation can reproduce the frequencies measured 
on a real toroid to a very high accuracy (Schliesser et al. ( 2008a[ )). Indeed, 
the average deviation between measured and simulated frequency is, on av- 
erage, below 2%. This emphasizes the reliability of the FEM. Note also that 
probing of the modal displacement patterns using a scanning probe technique 
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Fig. 10. Resonance frequencies of silica microspheres for the (n, /) = (1,0) (upper 
branch) and (n, /) = (1,2) (lower branch) spheroidal resonance modes ( |Ma et aT 
( |2007| )). Curves are expected resonance locations for free boundary conditions, fol- 
lowing an inverse dependence on the sphere radius with Vt^/2Ti ^ 2280m/s/i? for 
the / = and Vt^/2Ti ^ 1580 m/s/i? for the / = 2 mode. Symbols represent measured 
spheres. The panels on the right show the corresponding displacement patterns as 
obtained from the FEM. The left and right halves of the sphere show the displace- 
ment of the volume elements at two different times separated by half the oscillation 
period. 



has confirmed the shape expected from simulations in an earlier experiment 



dKippenberg et al.| ( |2006| )). 



Various mode families with qualitatively different displacement patterns are 
recognized. Some modes involve mainly motion of the silicon pillar (number 6, 



13 and 18 in figure 12) and therefore couple only negligibly to the optical mode, 
rendering it irrelevant for the purpose of cavity optomechanics. A torsional 
mode (number 4) is also observed, but it couples only weakly to the optical 
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modes (cf. section 3.3). Modes number 1, 3, 5, 7, 11 and 16 are characterized 
by an n-fold symmetry under rotation by an angle vr/n, where n = 1 ... 6 in 
this case. These modes, which we refer to as "crown" modes, involve mainly 
sinusoidal oscillation of the toroid in the z-direction. They typically follow a 
quadratic dispersion relation, that is, their eigenfrequency is proportional to 



the square of the number of nodes along the circumference (Schliesser et al 



( 2008aD ). Modes of this family are doubly frequency degenerate as for each 
mode in this family there exists a second mode in which the positions of nodes 
and antinodes are swapped. In measurements on real toroids, this degeneracy 
is often lifted due to small deviations from perfect rotational symmetry of the 
structure, leading to a splitting into a resonance doublet. 



In the context of cavity optomechanics, the most interesting modes coupling 
strongly to the optical field are the ones with (nearly) radially symmetric 



radius 




major radius 




Fig. 11. Geometric parameters required to describe a sphere (right) and a toroidal 
WGM resonator, where silica is shown in blue and silicon in grey. 






Fig. 12. Displacement patterns of the 19 lowest-frequency modes (see text) of a 
toroidal microresonator as calculated using the FEM, indicated both as the de- 
formed shape and in the color code (increasing displacement from dark blue to dark 
red). Mode number 4 involves mainly torsional motion, which is illustrated by plot- 
ting the displacement of originally parallel slices through the structure. Figure from 



ref. (Schliesser et al. (2008a)). 
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Fig. 13. Comparison of simulated and measured frequencies of a toroidal silica 
microresonator. The relative deviation is below 2% on average. Figure from ref. 
dSchliesser et al.| ( |2008a"| )). 



displacement patterns (number 2, 8, and 14). One usually distinguishes flex- 
ural modes, in which the displacement is mainly along the z-axis (number 2 
and 8), and the radial modes involving mainly radial displacement. Strongest 
coupling is caused by the fundamental radial mode, which we refer to as the 



radial breathing mode (RBM) with number 14 in figure 12 The RBM in the 
measured toroid has a frequency of 75.1 MHz, and usually lies between 30 and 
120 MHz for typical toroidal geometries. 



3.2.4 Mechanical dissipation 

Another very important property of the RBM is its damping rate. As in the 
case of the optical modes, diflFerent dissipation mechanisms lead to a release 
of mechanical energy stored in the RBM to other degrees of freedom. The 
mechanical damping rate quantifies the rate at which this takes place, and 
is again expressed in terms of a quality factor 

Qm = f^m/r^. (78) 

The different damping mechanisms all contribute to the total damping, and 
the resulting quality factor can be written 

Qm = Qgas + Qclamp + QtLS + Mother (^9) 

with contributions by the surrounding gas (Qgas)? by clamping losses {Qcilmp)^ 
by two-level systems (Qtls) by other damping mechanisms (QotLr)- 
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In the initial optomechanical experiments with silica microtoroids the mechan- 



ical quality factors reached values of 3,000 (Kippenberg et al. (2005); Schliesser 



et al. ( 2006D ). By replacing the surrounding nitrogen gas with helium (which 



has lower viscosity rj and molecular mass M), the quality factors rose be- 
yond 5,000, clearly indicating that the losses were dominated by gas damping. 



From th e known drag on a sphere oscillating in a viscous fluid (Landau and 
Lifshitz (1987)) one can expect a scaling Qgas oc {r]pM)~^^'^ in this regime. 



Reducing the pressure p of the ambient gas should therefore raise Qgas- In- 
deed, at pressures above lOmbar the quality factor was observed to increase 
as Qgas oc indicative of the viscous regime, while for lower pressures, in 



the molecular regime, Qgas oc was observed (Anetsberger et al. (2008)). 
At pressures below 1 mbar, gas damping becomes irrelevant (Qgas 3> Qm) and 
the quality factor is observed to converge towards a saturation value at lower 
pressures. Among different toroids, there are signiflcant variations in the range 
between 1,000 and 30,000 for this value. This is due to clamping losses, which 
strongly depend on the geometry of the sample as is detailed below. 

To study this effect of clamping losses more systematically, a toroidal resonator 
was underetched in several steps, to reduce the radius of the silicon pillar, and 
thereby increase the relative undercut, which we deflne as 

. _ silicon pillar radius . . 

relative undercut = 1 — . (80) 

silica toroid major radius 

The quality factors and resonance frequencies of the RBM of six toroidal 



resonators on the chip were measured for each etching step (flgure 14). Strik- 
ingly, while the resonance frequency reduces monotonically, the quality factors 
vary in a non-monotonic manner, with a distinct minimum for an undercut 
of approximately 0.7. This behavior was highly reproducible and observed for 
different independent samples. 

A series of simulations with increasing undercut reveals that the resonance 
frequencies of the RBM and a radially symmetric flexural mode are crossing 
each other for this undercut, as they have a different undercut dependence. 
However, the actual resonance frequencies of the structure exhibit an avoided 
crossing. At the same time, the modal shapes of the RBM and flexural mode 
hybridize in the crossing region. These two facts imply a normal mode coupling 
between the RBM and the flexural mode. 



To corroborate this conjecture experimentally, the undercut in a different sam- 
ple was again systematically increased. For each etching step, a highly sensitive 
measurement technique (polarization spectroscopy, cf. section ^ allowed the 
determination of the frequencies of both the RBM (J^r) and also the adjacent 
flexural (f^p) mode. Indeed, the measured frequencies and quality factors can 
be simultaneously reproduced using a simple model of two coupled harmonic 
oscillators, the intrinsic frequencies and and quality factors and 
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Fig. 14. Normal mode coupling between the RBM and a flexural mode. Left panels: 
Measured quality factors of the RBM of six samples on the same chip as a function 
of relative undercut, and the resonance frequencies of one of the samples (dots). 
Lines in the lower panel are results of simulations, and show not only the frequency 
of the RBM, but also of radially symmetric flexural modes. In the region where 
the frequencies are similar, for a relative undercut of ^ 0.7, the quality factors 
are strongly reduced. Right panels: Evidence of normal mode coupling between 
the RBM and a flexural mode in another sample. The frequencies fl± and quality 
factors Q± of the experimentally observed modes (dots) can be reproduced using a 
model of two coupled harmonic oscillators (red and blue lines). The frequencies and 
quality factors of the uncoupled modes, corresponding to a pure RBM and flexural 
mode are assumed to depend only linearly on the relative undercut (dashed lines). 
Adapted from ( [Anetsberger et al.| ( [2008^ ). 

Qy of which are linearly dependent on the undercut. The frequencies Q± and 
quality factors Q± of the new eigenmodes are given by (Anetsberger et al. 

poos] )) 



1 



\{2 




(81) 



where the intermode coupling g] 
data in figure 



14 



is an adjustable parameter (figure 14). The 
can be fit using S'im/STr = 14 MHz. This rather strong coupling 
is attributed to the asymmetry of the structure in the axial direction due to 
both the offset of the toroid from the disk ([Kippenberg et alT (2005)), and the 



fact that the silicon pillar supports the disk only from below. 



Summarizing the previous observations, we conclude that the admixture of 
a flexural displacement pattern to the RBM reduces the quality factor of 
the latter. This can be explained by noting that the flexural modes induce 
axial displacement also in the region where the silica disk is supported by 
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the silicon pillar. The periodic oscillation in axial direction launches acoustic 
waves into the pillar, through which the energy of the mode is dissipated, 
thereby deteriorating the quality factor. 

This qualitative understanding is quantitatively supported by an analytical 
model (Anetsberger et al. ( 2008D ), in which the acoustic energy loss is esti- 
mated as the power radiated by a membrane of area (area of the silicon 
pillar) oscillating with an axial displacement Az{r) at frequency 



mech 



2 ^2 



d r. 



(82) 



where v is the sound velocity and p the density. For geometry parameters close 
to the modal crossing, it was indeed found experimentally that 



Qclamp ^ 



m-^mech 
mech 



(83) 



where E'mech is the total mechanical energy of the mode. Advantageously, both 



^mech and Pmech cau bc simulated using FEM (Anetsberger et al. ( 2008[ )). As 
an aside we note that an expression similar to (82) and (83) was obtained in a 



rigorous theoretical analysis based on a phonon tunneling approach (Wilson 



Rae ( 2008D ). To reduce Pmech, two strategies are immediately evident: Either 
the clamping area is minimized, or the axial displacement Az{r) in this 
region is reduced. Minimizing the clamping area is possible, for example, by 
fabricating toroids with a very strong undercut, as shown in figure [T5| in which 
the silica disk is supported by a "needle pillar" of sub-micrometric diameter 



(Schliesser et al. (2008b)). 



Alternatively, by introducing spokes into the silica disk, it was possible to 
strongly reduce the coupling of the radial motion to an axial displacement in 
the clamping region (Anetsberger et al. ( 2008] )). Both presented approaches 
yield good results with quality factors exceeding 30,000 at frequencies around 
40 MHz dSchliesser eFal] ( |2008bD ; [Anetsberger et al.| ( |2008D ). 



In figure [16] we show an overview of mechanical quality factors achieved in 
typical samples with different major diameters, all strongly undercut (> 90%) 
and measured in vacuum. A clear trend to higher quality factors for larger 
cavity sizes (and therefore lower frequencies) is observed. The scatter in the 
data of neighboring toroids (with very similar refiow preform and final pillar 
shapes) indicates that clamping losses depend sensitively on geometry parame- 
ters. Finally, we note that we have consistently observed higher quality factors 



in disks than in toroids. This is attributed to a reduced offset (figure 11) of 



the oscillating mass, which, in toroids, mediates the coupling of the RBM to 
radially symmetric fiexural modes with strong dissipation to the pillar. 



The highest quality factors which were experimentally achieved at room tem- 
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Fig. 15. Silica toroidal resonators with ultralow mechanical damping, a) Reduced 
clamping loss by supporting the silica disk with a "needle" pillar. The central and 
right panels show the tip of an intentionally broken pillar, which has a sub-micron di- 
ameter. A quality factor of 30,000 is reached with such structures for a 40 MHz-RBM 



(Schliesser et al. ( 2008b| )). b) Reducing clamping loss by engineering the axial dis- 
placement in the clamping region using a "spokes" design. Such resonators achieved 
mechanical quality factors up to 32,000 at 38 MHz, and 50,000 at 24 MHz at room 
temperature ( [Anetsberger et aL ( 2008| )) limited by TLS losses. 



perature are Qm ^ 50,000 above 20 MHz. For these structures however, sim- 
ulations clearly indicate Qciamp ^ Qm according to equation (83), indicating 



that a different dissipation mechanism must now be dominant. A strong tem- 
perature dependence of the quality factor (allowing for values up to 80,000 



at 110°C) suggests a temperature-dependent dissipation mechanism (Anets- 
berger et al.] ( |2008D ). 



Indeed, the intrinsic damping of acoustic excitations in glass is known to follow 
a universal behavior observed in many amorphous solids (Pohl et al. ( 2002D ). 
This effect is attributed to the coupling of strain fields to structural defects 
in the material. While the microscopic nature of these defects is not precisely 



32 



40 Mm, 106 MHz 



O 



o 
E 





6000 
















4000 






ll 










T ir 2000 


' ^ i ' i 




s 





sample # 



80 |jm, 48 MHz 









20000 






















1 ouuu 
























10000 - 










.1 


T 





1 30 

sample # 



100 |jm, 42 MHz 









=1 

z 


15000 


I 






pi 0000 








5000 





1 26 

sample # 



100|jmdisk, 37.4 MHz 

T ^ 25000 

20000 ^ 

1 5000 
10000 
5000 



1 9 1 4 

sample # sample # 

Fig. 16. Overview of mechanical quality factors of the RBM measured in typical 
toroids. Major diameters of the pre-reflow disk and the approximate resonance fre- 
quency are indicated in the figure captions. Each bracket indicates the span within 
which the mechanical quality factor was found in a pair of measurements, in which 
the probing laser was red and blue detuned (dynamical backaction modifies the 
measured effective mechanical quality factor in opposite directions in the two cases). 
Groups of brackets belong to the same toroid. The last panel shows measurements 
taken on silica disks prior to the reflow. 



a 

o 
"c 

03 

o 

CD 



130 Mm, 30.5 MHz 



15000 
10000 
5000 3 



If 



known, it can be successfully modeled by a distribution of effective two-level 
systems (TLS) with two stable equilibria, represented by a particle in an asym- 



metric double- well potential (Jackie (1972)). These potentials are character- 
ized by the energy asymmetry A of the ground states in both potentials (the 
two levels involved) , and the height V of the energy barrier separating the two 



wells as shown in figure [T7| Oscillating strain fields associated with an acous- 
tic excitation modulate the energy asymmetry A between the two potential 
minima, and thereby couple to the TLS. 



To a very good approximation, the resulting Debye relaxation of acoustic 



excitations leads to a quality factor given by (Vacher et al. (2005)) 



Qtls — 



/+CX) r+oo / 
/ P(A, V) sech^ ; 



A 



2kBTj 1 + fi2r2 



dVdA, (84) 
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Fig. 17. Double-well potential used to model the structural defects in glass. The two 
individual wells (dashed lines) are usually assumed to be identical, but to have a 
ground-state energy that differs by an energy asymmetry A. The two stable equi- 
libria are separated by a barrier of height V. 

where 7 = ^dA/de is the change of the potential asymmetry as a function of 
strain 5, p is the density, v sound velocity, and P(A, V) is the distribution of 
TLS in the energy parameters A and V in the sense that P(A, V)dVdA is a 
volume density of defects. 

At room temperature, the relaxation between the two stable states is domi- 
nated by thermally activated processes. The relaxation time is thus given by 
an Arrhenius-type law 



Following the arguments given in reference (Vacher et al. ( 2005| )), a sensible 



choice of the distribution function P(A, V) yields eventually a quality factor 
of 

^--^ = ^-''^[-^)k^Jo [vj ' " l + QVoVW ^^ (86) 

reproducing the experimental data over four orders of magnitude in frequency 
(11 kHz . . . 200 MHz) and two orders of magnitude in temperature (from a few 
Kelvin to above room temperature) for the parameters Vq = (667 ± 21) K • /cb, 
C = 0.28 ± 0.03, logio(ro/s) = -12.2 ± 0.18 and Vq/A, = 7.7 ± 0.7. 

While usually measured as the attenuation of large-amplitude acoustic waves 
in bulk material, the very same temperature dependence of the quality factor 
was found for the REM of spokes toroids with sufficiently low clamping losses 
(Arcizet et al. ( 2009| )). Figure 18 shows the measured quality factor of the 



RBMs of two samples at frequencies of 36 and 63 MHz. Simultaneously with 
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the damping, relaxation of the TLS also leads to a change in the sound velocity, 
giving rise to a frequency shift of the mechanical modes (Vacher et al. ( 2005[ )), 



-^•erf 

2 V Ac / ksT Jo 



1 



1 + n^rie^v/k^T 



dV, 
(87) 

which is also shown in figure [TS} At temperatures above 10 K, the damping 



can be accurately described by equation (86), with a peak damping at about 



50 K leading to a minimum quality factor of Qtls ^ 500. 
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Fig. 18. Quality factor of the RBMs of two samples which are dominated by the 
damping due to coupling to structural defects in glass. Solid line is a fit using 



equation (86), and the dashed lines are experimental data from an acoustic wave 
attenuation experiment at 40 MHz (Bartell and Hunklinger ( 1982[ )). The inset shows 
the relative frequency shift of the RBM induced by the relaxation of TLS. From ref. 
dArcizet et aLl ( |2009D ). 



Below 10 K, the relaxation is dominated by tunneling processes between the 
two equilibria (Jackie (1972); Tielbiirger et al. ( 1992[ )), instead of the thermally- 
activated relaxation. The tunneling relaxation is responsible both for the low- 
temperature plateau (Qm ^ 1200 at 5 K) and the roll-off at very low tempera- 
tures with Qtls ^ For completeness, we note that apart from damping 
via relaxation, at sufficiently low temperatures, direct absorption of acous- 
tic waves by the TLS also leads to damping, which saturates at high enough 



amplitudes (Jackie (1972); Hunklinger et al. (1973)) 



Other damping mechanisms, as described by Qother, include thermoelastic 
damping (Zener (1937, 1938])), damping by anharmonicity (Vacher et al. ( 2005[ )), 
and surface effects (Ekinci and Roukes (2005)). For silica microtoroids or 
spheres, these effects are individually estimated to lead to limiting quality 
factors on the order of 10^ or more. In particular when operating the res- 
onators at cryogenic temperatures — as required for advanced experiments in 
cavity optomechanics — these damping mechanisms can be safely neglected 
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compared with the damping due to TLS. 
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3.3 Optomechanical coupling 



In order to describe optomechanical coupling in microspheres and microtoroids 
introduced in the previous section, it is advantageous to map the mechanical 
modes of interest to an effective one-dimensional mechanical oscillator, de- 
scribed by a displacement x, which parametrically modulates the optical reso- 
nance frequency through a non-zero = duc/dx. For optomechanical devices 
which host optical and mechanical modes with complex three-dimensional 
mode distributions such as silica microtoroidal resonators (figure 19), this 
mapping can be non-trivial. 

optical WGM mechanical RBM 





Fig. 19. Mode shapes of optical and mechanical modes in a silica microtoroidal 
resonator (FEM simulations), which have to be taken into account to quantify the 
strength of optomechanical interaction. 



3.3.1 Mapping to a scalar displacement 

To formally describe this mapping, a weighting function w{f) is introduced, 
mapping the displacement field u{f^ t) to a scalar displacement x according to 

x{t) = j w{r) • u{r,t)d^r = Y.^n{t) j w{r) • <(f)rfV (88) 
= J2c^{t){w,if\, (89) 



where the decomposition of the displacement pattern u as presented in eq. (61 ) 
is used. For each mechanical mode, the overlap integral {w^ vP^) determines the 
relative strength of the transduction of the n-th mode amplitude c^(t) into 
the optically sampled displacement x{t). The global normalization of w{r) is, 
in principle, arbitrary, as it can be compensated by an adequate choice of 
g^ = dujc/dx. For example, it can be chosen such that the displacement of 
a particular part of one mode's displacement pattern (such as the antinode 



of a beam's fundamental mode) directly corresponds to x (Eichenfield et al. 
( |2009bD ). 



For the most generic optomechanical systems, there are other obvious choices. 



37 



For example, for a Fabry-Perot cavity, the normalization is chosen such that 
X corresponds to the center-of-mass mirror movement if it was displaced as 
a whole. Thus, if the laser spot on a mirror at z = Zq, which contains the 
mechanical modes, is given by a rotationally symmetric Gaussian with a waist 



of W{)^ the weighting function reads (Pinard et al. 



(1999); Braginsky et al. 



w{r) = 



TTWq 



(90) 



With this weighting function, the physically correct frequency shift is obtained 
using the coupling constant = —uj^/L^ where L is the total length of the 
cavity. 



In the case of silica WGM resonators, with their three-dimensional distribu- 
tions of optical and mechanical fields, the calculation of the effective displace- 
ment is more difficult. A sensible approximation may be obtained by consid- 
ering the mechanical displacement as a perturbation, which does not modify 
the optical fields, but only displaces polarizable matter within the optical field 
distribution. The resulting relative frequency shift equals the relative change 
in the electromagnetic energy stored in the mode, as the number of stored 
photons, each of energy /kJc, is conserved in the cavity (Arnold et al. ( 2003D ). 
The resulting frequency shift is 

5uj, _ 5E,^ _ J lE\f)Vs{r) ■ u{f)dh 

uj, E,^ 2!\E\r)e{f^(Pr ' ^ ^ 

where £'^(r) is the squared magnitude of the electric field, and the factor 
2 in the denominator is due to the fact that equal amounts of energy are 
stored in both the electric and magnetic fields. Note however that the magnetic 
energy does not shift upon a displacement of the boundary, as the magnetic 
permeability of glass is very close to unity. For a resonator made out of a 
homogeneous dielectric material, the integral in the nominator is essentially a 
surface integral over the boundary, as Ve{r) is zero everywhere except for the 
surface. We use the coupling constant 5^0 = —^c/R^ expressing the frequency 
shift as a consequence of an effective radius change. One then finally obtains 

w{r) = . - E\f)Ve{f). (92) 



As equation (92) is difficult to analytically evaluate in complex geometries, a 



useful approximation is given by 

~ T^^i^ - ^o)5(r - R)er. (93) 

where R is the major radius and z = Zq the plane of the equator of the 
toroid. This weighting function essentially considers the transverse size of the 
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optical mode as negligibly small compared to the scales of the displacement 
patterns, and the resulting displacement x corresponds to the change of the 



cavity radius, which is sampled by the optical mode. In this work, equation (93) 



is applied to derive the displacement induced by the excitation of a particular 
WGM from the results of FEM simulations. 



We finally note that in a dielectric resonator, it is important to also consider 
strain-optical effects, that is, a strain-dependent refractive index leading to 
additional resonance frequency shifts for a given excitation of the mechanical 
mode. This effect was found to dominate the optomechanical coupling in a 
cryogenic sapphire microwave WGM resonator ([ Locke and Tobar| ( |2004D ). To 
assess the relative contribution in silica WGM microresonators, we may use the 



analytic expressions for the strain field (67)-(69) in a microsphere to calculate 



the corresponding change in the refractive index as experienced by the optical 
mode. Due to the homogeneity of the strain fields on the scale of the optical 
mode cross section, it can be well approximated by just evaluating the strain 
fields at the edge of the sphere using dllchenko et"aLl ( fl998l )) 



5 {n = P2Srr + Pisee + P2e4>4> 

8 (^~Otm ^ Pl^rr+P2S0e+P2S^^, 



(94) 
(95) 



where the required coefficients of the photo-elastic tensor are given by pi = 
0.121 and p2 0.270 ( Dixon] ( 1967D ). This leads to an extra frequency shift 
of about 30% (TE modes) and 50% (TM modes), as compared to the shift 
induced by the displacement of the boundary alone. Evaluation of the strain- 
optical coupling in a toroid is not possible analytically. For typical torus ge- 
ometries as used in this work, however, we can extract an extra frequency shift 
of less than 20% from FEM simulations, in agreement with an earlier estimate 
dKippenberg et al.| ( [20051 )). 



3.3.2 Effective mass 



As a global coupling coefficient go is used to quantify the coupling strength 
of an effective displacement x to the resonance frequency ujc of the cavity, it 
is necessary to absorb the different coupling strengths of different mechanical 
modes into another parameter. This parameter is referred to as the effective 



mass of the individual mechanical modes (Gillespie and Raab (1995); Pinard 
eFal] ( [19991 )). 



An operational definition of the effective mass of one particular mode (labeled 
in the following with an index n) can be derived from its potential energy Un 
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which can be recast from equation (76) to 



Un = ^Mr,nnMt)Y (96) 



using the free-boundary conditions (Pinard et al. ( 1999D ) and the definition of 
the moving mass 

M„= / pKif^l'd'r. (97) 
Jv 

As Cn{t) is experiment aUy not accessible, we want to express the potential 
energy in terms of the measured displacement of the mode, 

Xn{t) = w{f) . Un{r, t)dh = Cnit) (w , , (98) 

and therefore require 

Un = \m,s,^^l{xr,{t)f. (99) 
This immediately leads to the formal definition 

In practice, equation ([99]) is used to calculate the effective masses from exper- 
imental data (where, for a given ^'o, ^ is directly measured, and Un ^ k^T /2) 
or FEM simulations, from which Xn and Un can be simultaneously extracted. 

The eflFective masses of the fundamental sphere modes can be calculated an- 



alytically, as both the potential energy (77) and the resonance frequency 



(75) as a function of the radius R are known, yielding the numeric relation 
meff = 8470 kg/m^ -i?^ for silica, i. e. approximately 30 ng for a 30 /xm-diameter 
sphere. The numerical values for the RBMs in silica microtoroids are lower, in 
the range of 3 to 20 ng for the typical dimensions used. 



3.3.3 Forces acting on the mechanical modes 



To calculate the radiation pressure force acting on the mechanical modes, it 
is necessary to assess the momentum transfer from the optical mode to the 
dielectric medium. The fundamental starting point for such an analysis in 
a complex geometry such as a silica WGM resonator is the fiux of momen- 



tum density of the electromagnetic field in the medium (Landau and Lifshitz 
( [T984| ; [Pfeifer et al.| ( [2007| )) 



T = 



-EaDa — HaBa 





(101) 
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where H and B denote the usual electric and magnetic fields, indices i 

and j denote the cartesian components of vectors and tensors, and 5ij is the 
Kronecker-delta. In vacuum, the entity Tij is usually referred to as Maxwell's 
stress tensor. The body force /, that is, the force density acting on the medium, 
is given by the divergence of this fiux, plus a contribution from a temporal 
change of the fiux density 



BT 

f — _ \^ 



d_l_ 



Si 



(102) 



where 

S = ExH (103) 

is the real-valued Poynting vector. With D = eE and B = ixH we can write 
this as 



f = -^EEVe-^HHV^ + 



'^-l)l>^x#) (104) 

^SolLio J c^dt 



where MaxweU's equations V-D = 0, V-^ = 0, VxE* 



-/LidH/dt and 

V X H = +£dE/dt were used (Landau and Lifshitz (1984)). We will disregard 
in the following the last term, the so-called Abraham force ^ as it is usually 
small, and oscillates at the optical carrier frequency, so that it averages out on 
the timescale of a mechanical oscillation. Furthermore, we neglect the second 
term due to the close-to-unity magnetic permeability of normal glass and 
obtain finally 

/(f) -^E^ry V5(f). (105) 



This body force can be recast to 



E. 



f{f) = w{f)——- = w{f)F^ 



R 



rp 



(106) 



using equation (92) and the scalar radiation pressure force from equation (39) 



E. 



(107) 



Note that equation (106) also holds for the simplified assumptions leading 



to the simplified weighting function (93): If all optical power is concentrated 
to the rim of the toroid, an estimate for the force density can be derived 
from a simple consideration: A line element Rd^p contains the fraction dLp/2T{ 
of the total number of Eq^/Hu) intracavity photons, where E^^ is the total 
electromagnetic energy stored in the resonator. Within the time nRdip/c^ 



^ The nature of this force has also remained a contentious issue for decades (Pfeifer 
eFaD([2007|)). 
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each photon transfers a fraction of sin dip ^ dip of its momentum hk to the 
waU, yielding eventuaUy 

/(^) -^^^^^^ — -ro[z — ZQ)o[r — R)er w[r) (108) 

^ RdLp2^ huj nRdLp/c ^ ^ ^ ^ ^ ^ R ^ ^ 



Irrespective of the detailed form of the weighting function, the energy of the 
mechanical system is (Pinard et al. ( 1999D ) 



(109) 



where strain-optical effects are neglected for simplicity. This leads to the fol- 
lowing equations of motion for the mechanical mode amplitudes 



or, in the Fourier domain 

c„(r2) 



M^ini-n^-inVr.) 



it) 



(110) 



(111) 



where the damping r„ and the corresponding fluctuational force 5Ft have 
been introduced as well. As only the projection Xn of the excitation of the 
mechanical mode is measured, one obtains 



Xni^) = {w,U^)Cn{yt) 



((/(1]),<)+5Ft,„(Q 



1 



Frp(r2) + 



(112) 



The fluctuational thermal force SF^^n used here obeys (5FT,n(^) ^^T,n(^0) ^ 
M^, but the commonly employed approach is to introduce an eflFective thermal 
force SFth^n with replaced by meff,n leading to the form (55), and yielding 
finally 



the one-dimensional description used in section [2} 



(Frp(f])+5Fth,n(f^)) 



(113) 



As an illustration of the influence of the effective mass, flgure [20] shows the 
spectrum Sxx{^) of the fluctuations of x = for a toroid driven by 

thermal noise according to equation (113). Data were extracted from FEM 
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simulations for toroids the major radius of which has been continuously swept 
from 35 to 100 /xm. Frequencies are a direct simulation result, and the effective 



masses were extracted using equation (93). All quality factors were, for sim- 



plicity, assumed to equal 5000. Clearly, the strong signature of the RBM can 
be discerned for its low effective mass, on top of the background of weaker flex- 
ural modes. At a major diameter of ^ 75 /xm, the avoided crossing discussed 



in section [3. 2. 4| is also apparent. 

100 



CD 




50 60 70 80 90 
Major diameter (/im) 

Fig. 20. Color-coded displacement noise spectrum of a toroid with a 1 fim thick 
silica disk, a 4/xm minor diameter, 90% undercut and a varying major diameter, as 
simulated using FEM. The strong signature is from the RBM, the other traces are 
from flexural modes. 
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4 Ultrahigh-sensitivity interferometric motion transduction 



As a first application of optomechanical coupling present in WGM resonators, 
we consider monitoring of mechanical displacements using the optical degree 
of freedom. While in early work the separation of two resonators was measured 



in this manner (Ilchenko et al. (1994)), we focus here on the measurement of 
displacements related to the internal mechanical modes of WGM resonators 



trated in figure 21 



(Schliesser et al. (2008a)). The principal idea of such a measurement is illus- 



In the following, the theoretical limits, possible experimen- 



tal implementations, and experimental results will be presented. 




o 
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Fig. 21. High-sensitivity displacement sensing using optomechanical coupling. Left 
panel: an input field S[-a is launched into the taper and resonantly coupled to the 
WGM. The properties of the intracavity field a and the field 5out coupled back to 
the fiber taper are modified by the displacement x{t) of the RBM. Right panel: the 
launched field at frequency (jJ\ (green line) acquires a phase shift that depends on the 
mutual detuning of laser and cavity resonance frequency (lower red curve, dashed 
red curve indicates the WGM Lorentzian as a reference). If the displacement x{t) 
modulates the WGM resonance frequency, the phase (f{t) of the emerging field 5out 
is also modulated. This phase modulation can be detected with quantum-limited 
sensitivity by comparison with an optical phase reference. 



4-1 Theoretical limits of displacement sensing 



We first explore the theoretical limits in the sensitivity of the measurement, 
and restrict ourselves to the simple case of resonant probing A = (for a 
more general discussion, see section [7|. In this case, the dynamical properties 
of the mechanical oscillator are not affected by the presence of the light in 
the cavity, that is F^ba = ^dba = in equations (29) and (30). Still, the light 
can be used as a sensitive probe for the mechanical mode; in the following we 
discuss the performance and theoretical limitations of this method. 
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4.1.1 Quantum noise 



To analyze the fundamental quantum limits, we directly calculate from equa- 



tions (50)-(57) the noise in the light field at the output 



SsU^) = 5sU^) - .f^. (-zgoa5x{n) + + ^J^] (114) 



with an analogous equation for 55out(^)- For resonant probing, the mean field 
at the cavity output is real, and therefore the phase quadrature is directly 
given by 

Sqoutin) = I {-Ssoutin) + Ssl.iQ)) . (115) 

For the symmetrized noise spectral density^ of the phase quadrature at the 
output, defined by 

with 

27r(5(l] + n')S^f{n) = {SqU^)SqUn')) (117) 

we obtain 



with the correlation functions from equations (53)- (54). Evidently, the noise 
spectrum of the phase quadrature contains information on the mechanical 
displacement spectrum Sxx^ but also a background term (in this normalization 
equal to 1) which is due to the quantum noise. This background constitutes 
the fundamental imprecision of the measurement, and is given by 

if expressed as an equivalent displacement noise. In other words, the impreci- 
sion is the apparent displacement noise measured in such an experiment, due 
to the inevitable quantum noise in the measurement of the probing light's 
phase quadrature. 



Recast to experimentally more accessible parameters, equation (119) deter- 



mines the smallest possible displacement 5xmin(^) which can be measured 



Note that double-sided spectral densities are calculated in this chapter. 
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using a WGM resonator (Schliesser et al. (2008a)) 



femip(r2) 



A 



16% Trie 



Pin/^l^ 



1 + 



'or 



(120) 



where A is the wavelength in the medium and A/ the measurement band- 
width. Written this way, the importance of high-finesse cavities is directly 
evident: The smallest displacement that can be measured is roughly given by 
the wavelength, divided by the cavity finesse and the square root of the num- 
ber of photons accumulated in the measurement time (inverse bandwidth). 
The term 



'1 + AD? / is due to a less efficient transduction of the motion 
into phase shift for Fourier frequencies beyond the cavity cutoff (which can, in 



principle, be avoided using multiple cavity modes (Dobrindt and Kippenberg 
([20091))). 



As an aside, we remark that the same result can be obtained by considering 
the classical transduction of a displacement into the phase of the light exiting 
the cavity and comparing the result with the shot noise in the measurement 
process of the light phase ( |Schliesser et al. ( 2008b[ )). Importantly, equations 
(119) and (120) are independent of the particular strategy used to detect the 



light's phase, as long as it can be achieved in a quantum-limited manner. Two 
strategies to accomplish this, namely homodyne and polarization spectroscopy, 
are described in section 14.21 



As first discussed by Braginsky (Braginsky and Khalili (1992])), it is inevitable 
that the measurement of the oscillator's position disturbs it ("measurement 
backaction" ) . In the case of an optomechanical system, this is due to the 



fiuct nations of intracavity radiation pressure (Caves (1980j)), which can be 
written as 

5F,^{Q) = -hgoa (^Sa{n) + Sa^{n)) (121) 



in equation (50). Again from the known correlation functions, we obtain here 



Qba,qn 



Q2 + («;/2)2' 



(122) 



if the input noise is again only quantum noise. In this case, the force noise 



122 



122 



is referred to as quantum backaction. Evidently, equations (119) and 



fulfill the quantum-mechanically required inequality of the imprecision- 



backaction product (Braginsky and Khalili (1992)) 



— > — . 

4?7c ~ 4 



(123) 



By causing additional displacement fluctuations in the mechanical oscillator, 
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backaction noise also impedes the determination of the osciUator's displace- 
ment. The total uncertainty in the measurement is therefore given by 



where 



(124) 



(125) 



is the susceptibility of the mechanical oscillator. Obviously, a tradeoff in terms 
of the "strength" oc g^a^ of the measurement has to be made, as impreci- 



sion reduces, but backaction increases for "stronger" measurements (Caves 



(1980, 1981D ), as illustrated in figure 22 Optimum measurement conditions 
are reached for 



(126) 



or, equivalently an input power of Psql = ^'^i^O'Sql/^Vc- In this case, one 
obtains a total uncertainty of ( [Schliesser et al. ( |2008a} ) 

h 



(127) 



called the standard quantum limit (Braginsky and Khalili ( 1992[ ); [Caves ( 1980D ) 
in the case rjc = I. Its peak value is calculated at 



1 



h 



E 



(128) 



In this calculation we have explicitly considered the effect of the coupling con- 
ditions to the cavity, which can — as a unique feature — be varied continuously 
in the experiment by adjusting the gap between the coupling waveguide and 
the WGM resonator. The SQL is approached most closely in the overcoupled 
limit Tex <C tq. It is noteworthy that the fibre-taper coupling technique to 
microtoroids can deeply enter this regime, and 100 • Tex < tq {rjc = 99%) has 
been demonstrated ( |Spillane et al. ( 2003D ). 



4^ 1^2 Laser technical noise 



The previous derivation deals with the fundamental sensitivity limits. A fre- 
quent technical limitation is due to excess (beyond the fundamental) noise of 
the laser used for probing. Frequency noise in the laser, characterized by a 
power spectral density S;^^{Q)^ corresponds to a higher level of fiuctuations in 
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Fig. 22. Quantum limits in the measurement of mechanical displacements (blue 
line). For weak probing < ftgQL, measurement imprecision dominates the total 
uncertainty, while for stronger probing > ttgQL, the noise in the mechanical os- 
cillator induced by quantum backaction dominates the uncertainty. For optimum 
measurements with d^ = ^SQL' uncertainty is at the standard quantum limit 

S'x^^lf^m) = f^/T^eii^m^m' Under laboratory conditions, thermal noise is addition- 
ally present (red line). 



the input phase quadrature, 



Ms- 



|2 



^S(^) = l + ^^^-(^), (129) 



raising the background on top of which the displacement spectrum has to be 
observed. Note also that most schemes to measure the phase of the light rely 
on a phase reference. If this reference is noisy, because it is derived from the 
same noisy laser, the imprecision in the displacement measurement is given 

by 



ST/^'i^) = (130) 
if the frequency noise overwhelms quantum noise in the measurement. 
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4- 1.3 Thermorefractive noise 



Another important source of noise potentially preventing the measurement 
of mechanical displacements with quantum-limited sensitivity is fluctuation 



of the resonance frequency of the WGM ( |Gorodetsky and Grudinin (2004); 
Matsko et ah] ( |2007D ) which are not related to the mechanical oscillators. For 
a dielectric resonator as silica microspheres or -toroids, the dominant eflFect to 
be considered here are fluctuations of the refractive index due to temperature 
fluctuations. At any flnite mean temperature T, the actual average tempera- 
ture Tv in a volume V fluctuates according to (Landau and Lifshitz (1980)) 



CpVp' 



(131) 



where p is the material density, and Cp the speciflc heat capacity. This applies 
in particular also to the mode volume, within which the WGM samples the 
temperature-dependent refractive index. 

To calculate the frequency spectrum of the resulting fluctuations, Gorodetsky 



and Grudinin (Gorodetsky and Grudinin (2004)) have used a Langevin ap- 
proach, introducing fluctuational driving terms into the dynamic equations of 



temperature diflFusion (Braginsky et al. (1999, 2000)). As a result, the impreci- 
sion due to thermorefractive noise in a silica WGM resonator can be estimated 
to amount to 



R^Ssn/n{^) 




dq 



(132) 



where k is heat conductivity, R the cavity radius, and d and h the trans- 
verse mode dimensions. At low frequency 10 MHz), this noise indeed can 



dominate the measurement imprecision (cf. section 4.3) 



4-2 Experimental techniques 



4-2.1 Homodyne spectroscopy 

A commonly employed technique for quantum-limited phase measurement is a 
balanced homodyne receiver (Yuen and Chan (1983)), which has been used in 
earlier optomechanical experiments ( Hadjar et al.| ( |1999l ); Briant et al. (2003); 



Caniard et al. (2007a)). One possible adaptation of this technique to the ring 



topology of a WGM resonator is shown in flgure 23, The probing (or signal) 
beam and a phase reference beam, referred to as the local oscillator (LO) are 
derived from the same laser, in this case a monolithic Nd:YAG laser operating 



49 




(b) 



(c) 



Balanced 

homodyne 

detection 



Time (a.u.) 



frequency control 




Signal processing 



-50 50 
Detuning (MHz) 



Fig. 23. (a) Optical interferometric displacement transducer based on homodyne 
spectroscopy of light transmitted past the cavity ("/x-toroid"). The phase of the 
local oscillator is actively stabilized ("LO phase control"). Details are given in the 
text. PBS0-PBS3, polarizing beam splitters, (b) Cross section through the fiber 
taper and the toroidal rim in the coupling region. The polarization in the taper is 
slightly mismatched with the polarization of the cavity mode. Thus only the part 
E'cav of the total field couples to the WGM, the other component £'iock can be used 
for the stabilization of the local oscillator phase. The components ^'cav and £^iock cire 
separated in PBSl. (c) Signal in the balanced receiver for a scanning local oscillator 
(dotted, blue) at low power, and for the locked LO (red). The shown locked trace was 
recorded for about 5 seconds, (d) Typical experimental error signal in the balanced 
receiver when the laser is scanned over a cavity resonance with the local oscillator 



locked to the appropriate phase. Figure from ref. (Schliesser et al. (2008a)). 



at A = 1064 nm. This source exhibits quantum-limited amplitude and phase 
noise at Fourier frequencies Q/27r > 5 MHz and power levels Plo + P ^ 
5mW of interest. As tuning speed and range of this laser are limited, it was 
found advantageous to use a home-built external-cavity diode laser for pre- 
characterization of several samples until a suitable toroid was found. The 
probing beam is sent through the coupling taper and interacts with the WGM 
of the studied toroid. The LO travels in the reference arm of a Mach-Zehnder 
interferometer over the same distance. It is finally recombined with the signal 
beam at a polarizing beam splitter (PBSl). Spatial mode matching of the 
incident beams is enhanced by using single-mode fiber as mode filter on the 
local oscillator. After spatial recombination, interference is enforced using a 
retarder plate and another polarizing beam splitter (PBS2). 
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As the relative phase of the two interfering beams is also subject to drifts and 
fluctuations, due to, for example, temperature drift of the fiber in which the 
reference beam propagates, active stabilization is necessary here. In one pos- 
sible implementation, this is accomplished by purposely introducing a small 
polarization mismatch between the light in the taper region and the either 
predominantly TE- or TM-like WGM modes of the microcavity. The polariza- 
tion component of the probing beam which does not interact with the WGM 



resonance can then be used to stabilize the phase of the LO (see figure 23 ) . In 
this case, a dispersive signal, 

277c A 



^^^^ = A2;(^/2)2 V-PcavPLO, (133) 

is obtained at the output of the balanced receiver comparing the other polar- 
ization component (which is coupled to the WGM) with the LO. Here Pcav 
and Plo are the powers of the probing and local oscillator beams and h{A) 
is the power difference measured between the two employed receivers. Fig- 
ure 23 shows an example of an experimental trace obtained when scanning 
the laser detuning. Evidently, due to its dispersive shape, the d.c. -component 
of this signal can be used to lock the laser frequency to the WGM resonance 
frequency using electronic feedback. 

In this way, the mean detuning can be stabilized to A = and with A = 
A — go5x{Q) the signal is directly given by 



h^^,{dx, n) y i + (j^/(^/2))2' (134) 

where the reduced signal strength for Fourier frequencies Q beyond the cavity 
cutoff can be viewed as a consequence of the reduced buildup of intracavity 
sidebands ^ ^ . For a strong local oscillator Plo ^ ^cav, the detection noise 
is given by shot noise caused by the local oscillator beam, and the fiuctuations 
in the detected differential power are simply 5h ^ VPlo ^^i- Comparison with 



the signal (134) induced by displacements 5x then directly give the sensitivity 



derived in equation (119) 



Furthermore, as equation (133) evidently only depends on the mutual detuning 
of laser and WGM resonance, frequency fiuctuations of the laser are indistin- 
guishable from fiuctuations due to mechanical displacement. In the case of 
frequency noise of the laser, this leads directly to the imprecision described 



by equation (130). On the other hand, an intentional frequency modulation 
of controlled modulation depth of the laser can be utilized to calibrate the 
measured signals: a frequency modulation of Suj corresponds to a displace- 



A more detailed calculation is presented in the supplementary information of ref. 



(Schliesser et al. (2008b)). 
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ment of amplitude go5x^ independent of the detuning and coupling condi- 



tions ^ 


' ( 


Tittonen et al. 


(1999); 


Hadjar et al. 


(1999 


); 


Gorodetsky and Grudinin 


(2004 


); 


Schliesser et al. 


(2008b 


)). With the calibration at one particular mod- 



ulation frequency Omod, the measured spectra can be absolutely calibrated at 
all Fourier frequencies, taking into account the reduced sensitivity beyond the 
cavity cutoff at tv/2. 



4^2.2 Polarization spectroscopy (Hdnsch-Couillaud method) 

A simplified setup may be obtained by copropagating the local oscillator field 
in the same spatial, but orthogonal polarization mode as compared to the 
signal beam (Schliesser et al. (2008b])). Since the WGM modes have predom- 



inantly TE or TM character and are not degenerate, this guarantees that the 
local oscillator is not affected by the cavity. Due to common-mode rejection 
of many sources of noise in the relative phase between signal and LO (for ex- 
ample, frequency noise in the optical fiber), the passive stability is sufficiently 



enhanced to enable operation without active stabilization (figure 24). 



Enforcing interference between local oscillator and signal beams then corre- 
sponds to polarization analysis of the light (comprising both signal and LO) 
emerging from the cavity. While novel in the present context of a tapered fibre 
coupled microcavity, this is a well established technique to derive a dispersive 
error signal from the light refiected from a Fabry-Perot type reference cav- 



ity, named after their inventors Hansch and Couillaud (Hansch and Couillaud 
( [T9801 )). 



If fiber birefringence is adequately compensated, the error signal is 

2ri^K A 



^^^^ = A^TW^V^cavPLO, (135) 



identical to (133), and a typical trace is shown in figure [24](c). This is used 
to lock the laser at resonance A = with a bandwidth of about 10 kHz. 
Calibration of the spectra may be performed as described in the previous 
section. 



While this approach obviously allows one to reduce the complexity of the ex- 
periment, this arrangement proved more sensitive to slow temperature drifts 
in the polarization mode dispersion of the fibers employed, due to the large 
ratio of signal and LO powers, the magnitudes of which are only defined by 
the polarization state of the light in the fiber taper region. Improved stabil- 



^ We emphasize that for this relation to be valid for arbitrary modulation fre- 
quencies, it is necessary that the lengths of the two arms of the Mach-Zehnder 
interferometer are equal. 
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Fig. 24. Optical interferometric displacement transducer based on polarization spec- 
troscopy of light transmitted in the taper past the cavity ("/i-toroid"). (a) After 
phase modulation with an electro-optic modulator, the polarization is prepared 
with a first polarization control unit (PCUl). The cavity WGM defines signal and 
LO polarizations. A second polarization control unit (PCU2) compensates for fibre 
birefringence. Polarization analysis using a A/4 plate and a polarizing beam splitter 
enforces interference of the signal and LO fields, (b) Due to the polarization non-de- 
generacy of the WGM in the cavity, only one polarization component of the light 
interacts with the mode, (c) Typical error signal obtained when the laser is scanned 



over a cavity resonance. Figure from ref. (Schliesser et al. (2008a)) 



ity may be obtained by reducing the fiber length to its minimum of value 
of aproximately 0.5 m. For reasons of flexibility and convenience, the actual 
fiber length totaled to several meters in our experiment. Nonethel ess, sensi- 
tivities of 10~-^^m/\/Hz are achieved in toroids using this method (Schliesser 



et al. (2008b)). The intrinsic polarization selectivity of WGM renders the 



introduction of an additional polarizer, mandatory in the original implemen- 



tation (Hansch and Couillaud (1980)), obsolete. In an earlier experiment with 
a Fabry-Perot cavity (Hahtela et al. ( |2004[ )), the losses associated with an in- 
tracavity polarization element limited the finesse, and therefore the attained 
sensitivity to ^ lO"-^"^ m/v^Hz. 



4^2.3 Frequency modulation spectroscopy (Pound-Drever-Hall method) 



Another possible method to experimentally determine the detuning of laser 
and WGM resonance is frequency modulation spectroscopy as introduced by 
Pound, Drever and Hall (Drever et al. ( 1983D ) and discussed in great detail by 
Black ( pack] ( [2001] )). Fi gure 25 shows a possible application of this scheme 
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Fig. 25. Displacement measurement using the Pound-Drever-Hall method (Drever 
et al. ( 1983[ )). a) The phase of the probing laser is phase-modulated at a ra- 
dio-frequency of typically 50-100 MHz using a fiber-coupled electro-optic modu- 
lator (EOM). After interaction with the WGM, the detected signal is demodulated 
at the same frequency. The demodulation quadrature is chosen by adjusting the 
phase ip of the radio-frequency wave. In the simplest case, this can be accomplished 
by adjusting the length of the cable carrying the signal. The demodulated signal 
is low-pass filtered at a bandwidth well below the modulation signal. The resulting 
signal "PDH" can be used to monitor the detuning of laser and WGM resonance, 
and is also suited to electronically stabilize the laser frequency to the WGM. In 
addition, the transmission "T" of the WGM can be directly monitored. FPC, fiber 
polarization controller, b) Typical traces of transmission and PDH signals obtained 
with a silica toroidal WGM upon a laser frequency scan in the presence of thermal 
bistability. The satellite dips in the transmission signal appear when the modulation 
sidebands are scanned over the WGM frequency. 

to WGM resonators. The dispersive shape of the signal can be used to lock 
the laser to the WGM resonance using electronic feedback. Fluctuations of 
the Pound-Drever-Hall (PDH) signal beyond the feedback bandwidth then 
indicate fluctuations of the WGM resonance frequency with respect to the 
laser frequency, and can therefore be used to monitor displacements. 

While both homodyne and polarization spectroscopy can attain the funda- 



mental quantum-limited displacement sensitivity (119), the sensitivity of the 
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PDH method is reduced by a factor 1 + r^c + (1 — 27yc)^/2 where {3 is 



the phase modulation depth and Ji a Bessel function of the first kind ( Arcizet 



( [2007| )). Note that even for a maximally overcoupled cavity with Tj^ ^ 1 this 
yields 2 + l/2Jf In practice, the displacement sensitivity using this method 
is often limited by electronic noise in the detector. As essentially all light used 
in this scheme interacts with the WGM (there is no phase reference beam in 
a different spatial or polarization mode), the total power levels must be kept 
low in order to prevent strong thermal nonlinearities. These light levels (typi- 
cally a few microwatts) are not suflBcient to overwhelm the electronic noise of 
broadband light receivers. If available at the particular wavelength of interest, 
a low-noise optical amplifier such as an erbium-doped fiber amplifier (EDFA) 
can however ameliorate this drawback (Arcizet et al. ( 2009D ) at the expense 
of a higher noise figure. 



4.3 Observation and analysis of quantum and thermal noises 



Figure [26] shows data obtained using homodyne spectroscopy on a toroid of 
about 45 /xm radius. As long as the taper is retracted from the WGM evanes- 
cent field, quantum shot noise is observed to exceed the electronic detector 
noise. Note that while the detected shot noise (due to the local oscillator) is 
spectrally fiat (white noise) to a good approximation, the equivalent displace- 
ment noise plotted in figure 26 exhibits a calculated ^1 + Q'^/ (/^/2)^ frequency 
dependence beyond the cavity cutoflF at /^/2 27r • 17 MHz due to the Fourier 



frequency-dependent transduction (134) 



Approaching the taper to the resonator, and locking the laser to the WGM res- 



onance, a substantially diflFerent spectrum is observed (figure 26). The equiva- 
lent displacement noise is calibrated in absolute terms using an a priori known 



phase modulation at 36 MHz, as explained in subsection |4.2.1[ taking also the 
cavity cutoflF into account. While the background due to quantum measure- 
ment imprecision is at a level of 10~^^m/\/Hz at low Fourier frequency, a 
significantly higher equivalent displacement noise level is observed when cou- 
pling the laser to the WGM. 

The broad background particularly strong at low frequency can be quantita- 
tively reproduced by the model for thermorefractive noise (equation (132])), 
when no parameters except b and the absolute magnitude are adjusted by fac- 
tors of order 2. This is justified considering the approximations made in the 
derivation, and the incomplete knowledge on the transverse mode shape of 
the WGM probed in this experiment. We note here as an aside that thermore- 
fractive noise, measured here for the first time in a toroidal microresonator, is 
an important limitation for the generation of Kerr squeezing in these devices 



dRehbein et al] ( [20051 )). 
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Fig. 26. Equivalent displacement noise measured in a ^ 90 /xm-diameter silica 
toroidal cavity. Red, measured trace with laser coupled to a cavity resonance, in- 
cluding a peak at 36 MHz due to phase modulation for calibration purposes. Gray, 
measured shot noise with taper retracted from the cavity and black, detector noise. 
Models for mechanical noise (orange line) and thermorefractive noise (green line), 
and sum of the two models plus the shot noise background (blue line) are also 
shown. Figure from ref. (Schliesser et al. ( 2008a| )). 



On top of these backgrounds, a sparse spectrum of peaks is observed, which 
are due to different mechanical modes in the spectrum. In this measurement 
at room temperature, the thermal Langevin force largely dominates over ra- 
diation pressure force fluctuations. Therefore, each individual mode is driven 



by a random thermal force according to equation (113), and the thermal dis- 



placement noise spectra of the individual modes add up to the total measured 
equivalent displacement noise 

n 

where the symmetrized spectrum of the Langevin force is given, from equation 



(55), by 



S^pf{n) = hmee£,nTn^ coth 



2meff kBT, 



(137) 



where the second relation is valid as long as /cbT ^ HQ for the frequencies of 
interest. 



Figure 27 shows another example of a highly sensitive measurement, using 
the Hansch-Couillaud technique in this case. Beyond the clear signatures of 
the RBM at around 73 MHz, flfteen other peaks related to mechanical modes 
are observed. Zooming in on the individual peaks reveals that some are split. 
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Fig. 27. Broadband displacement noise spectrum recorded using the 
Hansch-Couillaud technique (top panel). Zooming in on the individual peaks 
(lower panels) reveals the precise frequency and linewidth of the modes (scale bar 
corresponds to a 100 kHz frequency span), some of which are split due to a lifted 
degeneracy. All observed peaks could be attributed to mechanical modes using 



finite element modeling (section 3.2.2). Figure from ref. (Schliesser et al. (2008a)). 



typically in modes whose degeneracy is lifted by residual asymmetry of the 
sample. The mode frequencies can be reproduced very accurately using finite 
element modeling; the peaks shown in this figure correspond to the modes 



discussed in section 3.2.2| (cf. also figures 12 and 13). 



Finally, in figure [28| we zoom in on a frequency interval that has the signature 
of a radial breathing mode (RBM) of a larger sample {R = 38 /xm). The 
measurement achieves a signal-to-background ratio of nearly 60 dB determined 
by measurement imprecision due to detection shot noise. This dynamic range 
exceeds the ratio 

We can therefore conclude that the imprecision background, at a level of 
1.1 am/ a/Hz, is below the standard quantum limit, calculated to \/s^^^{^m) = 
2.2am/\/Hz for this sample with n^/27r = 40.6 MHz, T^/27r = 1.3 kHz and 
meff = 10 ng. We emphasize however, that this does not imply that mea- 
surements with a better total uncertainty than the standard quantum limit 
are possible. Quantum backaction-induced fluctuations in the mechanical dis- 
placement increase the position uncertainty, but are masked by thermal noise 
in this measurement. 
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Fig. 28. High sensitivity- measurement of the RBM of a larger sample using the 
Hansch-Couillaud technique. Orange circles represent measured data, revealing also 
the signatures of two neighboring modes. The red line is a Lorentzian fit. Blue 
circles are the recorded measurement imprecision due to detection shot noise. An 
excellent signal-to-background ratio of nearly 60 dB is attained, corresponding to a 
measurement imprecision of 1.1 am/\/Hz, which is well below the standard quantum 
limit at 2.2am/\/Hz . 
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5 Observation of dynamical backaction 



In contrast to the resonant probing scenario discussed in the previous section, 
the dynamics of the mechanical osciUator is modified by radiation-pressure 
backaction if the optical mode is pumped in a detuned manner. This so-called 
"dynamical backaction" not only modifies the eflFective damping and spring 
constant of the mechanical oscillators, but, as we will show in the following, 
also leads to an energy exchange between optical and mechanical modes. 



Predicted as early as the 1960s by Braginsky (Braginskii and Manukin ( 1967D ), 
dynamical backaction has been observed early on in mechanical devices cou- 



pled to microwave resonators (Braginskii et al. 



Cuthbertson et al. 



( 1970 ) ; Braginsky and Manukin 
i )T996| ) ; |Locke et al.| ( |1998[ )). In 

the optical domain, dynamical backaction induced by radiation pressure has 
been first observed by the Vahala group in 2005 in the form of an oscillatory 



(1977 


); 


Blair et al. 


(1995 


); 



instability and studied in great detail 


((Kippenberg et al. 


( 


2005); Carmon et al. 


(2005); 


Rokhsari et al. (2005, 2006); 


Hossein-Zadeh et al. 


( 


2006)). 



Here, we present a systematic study of dynamical backaction as observed in 
silica microresonators, in particular its dependence on the relevant frequencies 
and A, and briefiy introduce the oscillatory instability described above. 
In the second part, we focus on the case of negative detuning A < 0. In this 
case, the light field extracts energy from the mechanical mode, leading to the 
cooling of the latter. This eflFect was first reported by our group at the MPQ in 
Garching (Schliesser et al. (2006)) and groups in Paris (Arcizet et al. (2006a)) 
and Vienna (Gigan et al. (2006)). Finally, we rule out thermal nonlinearities 
as the origin of optomechanical interactions in silica microtoroids. 



5. 1 Optical spring and optical damping 



For a detuned optical pump, we have found in section [2. 1.3| that the presence 
of light modifies the dynamics of the mechanical degree of freedom when it 
responds to an external force. In an intuitive picture, this can be understood 
as the consequence of the in-phase and quadrature response of the radiation- 
pressure force, when the mechanical oscillator is driven by the external force. 
The same result is formally attained using the quantum Langevin approach. 
Disregarding, in a first step, the quantum fiuctuations of the light {5s\^ = 
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<^Svac ^ 0) ^ , the radiation pressure force fluctuations in Fourier space are 



1 



1 



-i{A + n) + k/2 +i{A-n) + 



(139) 



equivalent to equation (25). As a consequence, the mechanical oscillator reacts 



to the thermal Langevin force with the effective susceptibility already derived 



in equation (27). 



To confirm these predictions, a series of measurements was taken using the 
setup described in figure [29j A 980 nm- wavelength diode laser was locked to 
a resonance of the silica microtoroid, simply using the transmission signal 
as an error signal, from which an oflFset can be subtracted to control the 
detuning. Moderate optical quality factors {Q < 10^) and low optical powers 
(Pin ^ 200 /xW) ensure that thermal nonlinearities are weak enough to still 
allow stable locking. 

Applying this procedure to both the red (A < 0) and blue (A > 0) wing of the 
optical resonance by changing the sign of the error signal, a detuning series 
can be recorded. From the transmission signal level with the laser locked to the 
side of the fringe, the relative detuning A/k. can be determined. At the same 
time, the fiuctuations of the transmitted power, as recorded by the spectrum 
analyzer, refiect the position fiuctuations of the mechanical modes. 

Driven predominantly by the thermal Langevin force with its essentially frequency- 



independent spectrum (137), the measured displacement spectrum directly re- 



veals the effective susceptibility of the mechanical mode. It is therefore possible 
to extract the effective damping and resonance frequency of the mode using 



the fit model (27). Figure 30 shows the data obtained from the 56.5 MHz- 



RBM of a silica microtoroid together with fits by the models (32) and (31). 



The measured changes in both damping and resonance frequency agree well 
with expectation. 



The resonance frequency shift is often referred to as "optical spring" effect 
(Sheard et al. ( 2004] )), as it originates from an optical restoring force propor- 
tional to the displacement of the resonator. It is interesting to note that this 
optical force can even exceed the natural restoring force of the mechanical 
oscillator, and thereby totally dominate the mechanical resonance frequency 



(Corbitt et al. (2007a)). For silica microresonators, due to the stiffness of the 



structure, this is typically not the case. In these devices however, another in- 
teresting effect occurs for narrow optical resonances (/^ < ^rn)'- in this case, the 
optical spring force changes its sign for small detunings, turning a restoring 



into an anti-restoring force and vice versa (Schliesser et al. (2006)), an effect 



^ This simplification is justified as long as the thermal Langevin force largely ex- 
ceeds force fluctuations due to the quantum nature of the light. 
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Fig. 29. Schematic illustration of the setup used for the measurement of dynamical 
backaction. A diode laser is locked to the side of the optical fringe by applying an 
electronic, computer-controlled offset to the transmission signal. This differential 
signal is used as an error signal in a feedback loop actuating the laser frequency 
by regulating both the position of the grating in the laser cavity, and the current 
pumping the laser diode. Once the laser is locked, the detuning is varied in small 
steps by adjusting the electronic offset. For each detuning, a trace from the oscillo- 
scope (to determine the actual residual transmission) and a noise spectrum from the 
electronic spectrum analyzer are taken. FPC stands for fiber polarization controller. 



not observed in other systematic studies of dynamical backaction (Arcizet 
et al.| ( |2006a| ); |Hossein-Zadeh and Vahala| ( |2007D ). 



The optically induced damping can provide both positive and negative damp- 
ing. For positive detuning A > 0, the total damping can reach zero. In this 
case, the mechanical mode, initially driven by thermal noise, starts to oscillate 
regeneratively. Specifically, solving Feff = for a threshold power, one obtains 



thresh 



X 



k/2 



k/2 



,(A - n^f + («:/2)2 (A + ^l^Y + {kI2)\ 



(140) 



for this optically driven mechanical oscillation to occur. This effect, often re- 
ferred to as parametric oscillatory instability (POI), has been reported for 
various systems, including silica microspheres (Ma et al. ( 2007[ )) with me- 
chanical modes at up to GHz- frequencies ( Carmon and Vahala| ( 2007| )). For 
light powers largely exceeding the threshold, nonlinearities neglected in the 
linearized models presented in this work lead to complex behavior such as 
multistability (Marquardt et al. ( 2006D ) and chaos (Carmon et al. ( 2007D ). 
For an in-depth theoretical discussion of the oscillatory instability, including 



also quantum effects, we refer the reader to references (Vahala (2008); Ludwig 



etgl] ( |2008a| ^ 



As an aside we note that in the regime n <C f^m = A discussed in greater 
detail in section [6} one finds the interesting relation 

i^thresh = 4y^PsQL. (141) 

This universally relevant power scale (for both dynamical and quantum back- 
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Fig. 30. Linewidth (a-d) and mechanical resonance frequency shift (e-h) of the me- 
chanical mode as a function of relative laser detuning A/hc for optical resonance 
linewidths n/27T of about 207MHz (a,e), 127MHz (b,f), 79 MHz (c,g) and 42 MHz 
(d,h). Dotted lines indicate the corresponding optical resonance Lorentzian, over 
which the laser was tuned. Dashed blue lines marks a detuning equal to the me- 
chanical oscillator's frequency of 56.5 MHz and full lines are fits from the models 
for dynamical backaction. 



action) is at the level of 30 fiW for typical parameters of silica microtoroids. 



5.2 Radiation pressure cooling by dynamical backaction 



In the preceding sections, we have only discussed the damping and resonance 
frequency of the mechanical mode and its modification by dynamical backac- 
tion. However, a major feature of light-induced damping is that it also changes 
the temperature of the mechanical mode. To introduce the concept of a "mode 
temperature" , let us first evaluate the amplitude of the displacement of a spe- 
cific mode, which is driven by the thermal Langevin force, by integrating its 
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noise spectrum over all Fourier frequencies 



r-\-oc _ rICi r-\-oo ^ _ , 

{Sx') = / SU^)^ = / \xm\'s%{n) 

J —oo ZTC J —oo 



•oo 

•CO 



27r 

2meffr^fcBT 



2^ 



-co ml^m^-^lf + mi)27:- 
This integral can be evaluated using the residue theorem, one obtains 



(142) 



(143) 



We may turn this result around and use it to introduce the mechanical mode 
temperature 

T^ = m,A{8x^)lkB. (144) 
With this definition, if the mechanical mode is only driven by the Langevin 
force, it ends up in thermal equilibrium with its environment, and = T. As 
an example, the root-mean-square (RMS) displacement {Sx^'^ of the RBM of 
a silica microtoroid is typically a few tens of femtometers at room temperature 



(figure 31). 




Time (|js) 



Fig. 31. Real-time recording of the displacement of a 40.6-MHz RBM of a silica 
microtoroid measured using the Hansch-Couillaud technique. A 2 MHz- wide band 
was filtered out of the displacement signal and selectively amplified to record the 
mechanical trace (red line). The gray line is a background trace recorded with the 
taper retracted from the cavity. Figure from ref. (Schliesser et al. ( 2008b[ )). 



In the presence of detuned pumping, the mechanical susceptibility is modified 
due to dynamical backaction, the thermal Langevin force, however, is not. If, 
therefore, the RMS displacement is calculated from the modified spectrum 

SU^) = We^m^^SfA^) (145) 

one obtains 

\mA{S^') ^ ll^k^T (146) 

Z ^ i eff 

as long as the mechanical oscillator can be described with its frequency- 
independent eflFective damping Fgff and resonance frequency Vt^^ for not too 
strong backaction eflFects . The mode temperature of the mechanical oscilla- 



See e. g. (Genes et al. ( 2008| )) for more general calculations 
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tor therefore is changed to 

T^ = ^T= T. (147) 

J- eff -1- m + -I- dba 

As for Fdba > one has < T the damping rate F^ba induced by dynamical 
backaction is often referred to as the laser cooling rate. 



Indeed it can be shown that F^ba is the rate with which energy is transferred 
from the mechanical resonator to the optical field. Returning, for simplicity, 
to the example of a sinusoidally oscillating boundary considered in section 



2.1| = Xo sin(r^ni^)), the cycle-averaged work done by the mechanical 



oscillator on the optical field can be calculated to 

2jr r^TT/Qm _ _ /I 



i^dba = F,^{t)x{t)dt ^ Fdba (^-mefff^^x^j (148) 

using only the elementary relations ([s]) and (11). 



In very general terms, this cooling effect arises by coupling the mechanical 
oscillator not only to the reservoir — consisting of all other mechanical modes 
present in the device, the gas and thermal radiation field surrounding it etc. — 
at room temperature, but also to the cooling laser field. Laser cooling therefore 
disequilibrates the oscillator with the reservoir, and brings it into a new equi- 
librium in which it is coupled both to the reservoir and the laser field. This 
field possesses an effective temperature very close to zero (see section [g] for 
the limitations), and acts as a "cold damper", by introducing dissipation, but 
only very little fiuctuations to the mechanical mode. 



We note here that the application of cold damping schemes has a long history 
in physics, and has been successfully applied in systems as diverse as electrom- 



eters and particle storage rings ( [Milatz et al. (1953); Milatz and van Zolingen 



(1953); Mohl et al. (1980)) (here often referred to as "stochastic cooling"^ 
Interestingly, cold damping has also been used to cool a mechanical mode of a 
mirror in a pioneering 1999 experiment at Laboratoire Kastler Brossel in Paris 



(Cohadon et al. (1999|)), and subsequently in many other experiments ( 


Arcizet 


et al. (2006b) 


; Kleckner and Bouwmeester| (|2006 ) 


; Weld and Kapitulnik 


(2006); 


Poggio et al. 


(2007);|Corbitt et al.|(|2007bD). However, these experiments all in- 



volve a complex hybrid electronic/optical feedback loop, whereas the method 
presented here relies solely on the intrinsic dynamics of radiation pressure. 



In figure 32 we show cooling results obtained on the RBM of silica microres- 
onators at frequencies around 57 MHz. In an experimental setup essentially 
identical to the one described in figure [29} a 980-nm wavelength diode laser 
was locked to the red wing of an optical resonance. Note that as the thermal 
bistability renders the red wing dynamically unstable under laser or cavity 
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frequency fluctuations (Carmon et al. (2004])), special care has to be taken in 
the implementation of the feedback loop stabilizing the laser frequency to a 
given detuning. 




Fig. 32. Cooling of the RBM of silica microtoroids. (a) Noise spectra of the RBM 
when a 50 MHz- wide WGM resonance is pumped at A ^ —Hi/ 2 with increasing 
power (0.25 mW, 0.75, 1.25 mW and 1.75 mW). The extracted increased damp- 
ing rates Feff correspond to mode temperatures Tm given in the legend, (b)-(e) 
Dependence of mechanical mode properties on the launched laser power (mea- 
sured on the RBM of a different toroid). The linewidth Feff and resonance fre- 
quency fieff vary linearly with input power in this range, and both the noise power 
Sxx{^)d^/'^^ ^ Tui and the resonance amplitude \/Sxx{^m.) vary linearly with 
the inverse linewidth F^ as expected. 

As the launched laser power is increased, the total damping increases, and 
correspondingly so does the width of the Lorentzian resonance. At the same 
time, the temperature of the mode is reduced by the optical pumping. 
Panel (a) shows four traces for detuned pumping of a WGM mode with A ^ 
—k/2 ^ —25 MHz. For the highest pump powers, the damping is increased 
beyond 450 kHz, and the corresponding reduced mode temperature is 11 K. 
The other panels show a systematic power series, in which a diflFerent torus 
was pumped at A —0.7/^ with powers between 20 and 200 /xW. 



5.3 Radiation pressure versus thermal effects 



In many early experiments, the optomechanical interaction has been mediated 
by thermal effects, sometimes referred to as "photothermal pressure" instead 



of radiation pressure (Hohberger Metzger and Karrai (2004); Gigan et al. 
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( |2006D ; [Harris et al.| ( [2007| ); |Favero etal] ( |2007D ; [Ludwig et al.| ( |2008bD ). In the 

following we present evidence for the fact that the optomechanical interaction 
in silica microtoroids is strongly dominated by radiation pressure due to the 
high (J^m/27r > 50 MHz) frequencies of the RBM. 



Prior work on the radiation-pressure-induced parametric oscillatory instabil- 
ity has provided independent evidence that radiation-pressure dominates the 



interaction between optical and mechanical modes ( |Kippenberg et al.| ([2005|); 
Rokhsari et al. ( 2006D ) in this system. For instance, it has been shown that the 
mechanical gain depends on the cavity finesse. This observation demonstrates 
that radiation pressure is responsible, since thermal forces only depend on the 
absorbed power and not on the cavity enhanced circulating power. 



5.3.1 Theoretical estimate 

To theoretically estimate the forces due to light absorption it is necessary 
to study the coupling of mechanical and thermal waves (or modes) of the 
structure. Thermodynamical considerations allow us to derive the mutually 
coupled differential equations for temperature and displacement distributions 



(Nowacki (1975)), adding a coupling term proportional to the linear expansion 



coefficient a to the heat diffusion equation and the equations of motion of the 



displacement ffeld (58). 



For small perturbations, we may however assume that the solutions of the 
uncoupled system {a 0) are still approximately valid, and heat transport is 
still dominated by diffusion. In this case, the resulting temperature gradients 
give rise to a thermoelastic body force ( [Nowacki] ( |1975D ) 

h,[r) = -(3A + 2/i)a Wr(f) (149) 

driving the mechanical modes (A and /x are the Lame constants). Importantly, 
this body force depends on the gradient of the temperature distribution. 

The effective scalar thermal force on a mechanical mode with a displacement 
pattern is determined by an overlap integral (/tet?^). Due to the very high 
mechanical resonance frequencies, the diffusion length for the temperature 



distribution Ad = y2k/cppQ becomes very short, for example about 50 nm 
for Q/27T = 50 MHz. If absorption takes place in the silica, we can therefore 
conclude that the temperature distribution is essentially given by the energy 
distribution of the optical mode, and the resulting temperature modulation 
in the volume can be estimated at 5T{Q) ^ 2Pabs/^Cpp]/mode5 typically < 
10 K/W for a typical 30 //m-radius toroid — a value confirmed by finite element 
modeling of heat diffusion in such a geometry. Due to the nearly symmetric 
temperature distribution in the radial direction with respect to the center of 



66 



the optical mode, the overlap integral over the gradient of the temperature 
distribution can be expected to yield only a small total contribution. For a 
rough estimate, we may use 



(AeO OC / VST{f)if^{f)d^r^ \if^{R)\27rRd^ dST/drdr 

JV JR-dm 

^ \if^{R)\27rRd^5T{R), (150) 

where rfm <C i? is the transverse diameter of the optical mode. Using 5T{R) / Pabs 
2K/W for the parameters described above (again confirmed by FEM), we ob- 
tain {fj.^il!^) / (ftQU^) ^ (9(10^). We note however that this result depends 
on the exact location of the heat source (absorption may also take place in 
a water or helium surface layer), and the cooling mechanisms provided by a 



surrounding medium, effects presently investigated in our group (Zhou et al 



( 2009D ). We also note that at cryogenic temperatures, thermoelastic coupling 



is weaker as the expansion coefficient drops below 2-10 ^atl.GK ( [White 



(1975); Arcizet et al. (2009)). 



5.3.2 Response measurements 



Experimentally, we have made response measurements (Rokhsari et al. (2005)) 
to quantify the different nonlinearities — due to the thermal, Kerr and radiation- 
pressure effects — encountered in silica microtoroids. To that end, two lasers 
at different wavelengths (980 and 1550 nm in this case) are coupled to two 
WGM resonances of a single toroid. One laser, referred to as the "pump", is 
amplitude modulated at a variable frequency while the other "probe" laser 
is used to measured the response of the WGM frequency to the pump laser 
power modulation. In the simplest case, this is accomplished by tuning the 
probe laser to the wing of a WGM resonance, and measuring the variation 
of its transmission at the same frequency f], most conveniently implemented 
using a network analyzer. Care is taken to suppress direct optical or electronic 
cross-talk of the pump modulation into this signal. 



Figure |33] shows the result of such a measurement on a 29 /xm-major radius 
toroid. Clearly, at low frequencies {Q/27r < IMHz), a strong modulation of 
the probe WGM frequency is apparent. This is due to absorption of pump 
light and the consequent modulation of the temperature-dependent expansion 
and refractive index of the toroid material. Above this frequency, a plateau is 
observed in the response, due to the modulation of the refractive index seen 
by the probe WGM, which varies with pump power due to the non-linear 
refractive index of silica (Kerr effect). Finally, around a mechanical resonance, 
a dispersive peak is observed. This is due to the excitation of mechanical 
modes by the modulated pump power. 
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To fit the response, we use the model 



with 



5cjprobe(^) = SuJthi^) + Su)k{^) + Su,p{Q) (151) 



\ n dr / \ 1 + zJz/JZi 1+ZJZ/JZ2/ C Tabs 



(152) 

5c.K(f^) = -c.c^^^^ (153) 

' '''"^fio(n). (154) 



For the thermal effect, we extract cutoff frequencies f^i ?^ 27r • 900 Hz and 
0^2 ^ 27r-69kHz. Furthermore, using the material parameters of silica and R = 
29 jam as well as ylgff ^2.5 /xm^, we can normalize all results to the measured 
Kerr response. In this manner, we extract /3i/rabs ^ 1-8 • lO^K/W/lOOns and 
P2/Tahs ^ 570K/W/100ns. All these parameters are very well in the range 
expected for the thermal effects. 



The combination of the radiation-pressure and Kerr responses gives rise to 
the dispersive signature around the mechanical resonance frequency: At mod- 
ulation frequencies slightly above the mechanical resonance, the mechanical 
degree of freedom oscillates out of phase and therefore counteracts the Kerr 
effect, which always reacts instantaneously to power changes. Furthermore, at 
the resonance frequency of the RBM we expect 



(^CJrp(f^m) 



'^^'"^^^ - 240 (155) 



with Fm = 15.7 kHz derived from the fit and the numerically determined effec- 
tive mass meff ^ 15 ng, in very good agreement to the measured value of 260. 
The extrapolated thermal effect, in contrast, drops more than four orders of 



magnitude below the observed mechanical displacement at (cf. figure 33). 
At mechanically non-resonant frequencies, this modulation would be largely 
dominated by a thermorefr active effect as compared to thermoelastic effects, 
as (dn/dT) /na > 10. Even if enhanced by the mechanical resonance, thermoe- 
lastic contributions may therefore be estimated at or below the 10~^-level. 



We finally note that thermally induced forces that may be related to the iden- 
tified thermal effects would be out of phase (by nearly 7r/2) with the driving 
pump laser modulation, since the relevant radio frequencies are well beyond 
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the thermal cutoff frequencies. As a direct consequence, the observed interfer- 
ence effects between mechanical displacement and the Kerr-induced modula- 
tion would form a single symmetric resonant peak, fundamentally different to 
the observed dispersive shape. This again confirms the dominance of radiation 
pressure in this system. 




10^ 10^ 10"^ 10^ 10^ lO'' 10^ 50 55 60 65 



Frequency (Hz) Frequency (MHz) 

Fig. 33. Pump-probe type response measurements of the nonlinearities present in 
toroidal silica microcavities. (a) Measured response (red dots) and fit (black line) 
modeling a two-pole low-pass thermal nonlinearity and an instantaneous Kerr effect, 
which becomes significant above 1 MHz. The data furthermore show a detector and 
cavity-induced cutoff above 100 MHz as well as an indication of mechanical reso- 
nances at 5.6 and 58 MHz, though not resolved. The green dashed line shows the 
extrapolation of the thermal effects to higher frequencies and below the measure- 
ment background at about 0.1 (gray dots). All data are normalized to the fitted Kerr 
response, (b) Higher-resolution data (red points) as recorded in the vicinity of the 
radial-breathing mechanical resonance at ^ 58 MHz. Full black line, model for the 
harmonic response plus a constant background due to the Kerr effect. Black and blue 
dashed lines illustrate the individual contributions from the harmonic response and 
the Kerr effect, respectively. Green dashed line indicates the extrapolated thermal 
response, and the green dotted line the expected response if the thermally induced 
displacement was resonantly enhanced. All data are normalized to Kerr response. 



69 



6 Resolved-sideband cooling 



The successful demonstration of cooling by dynamical backaction immediately 
raises the question of how strongly the temperature of the mechanical oscil- 
lator can be reduced. A series of theoretical and experimental investigations 
revealed the fundamental ( Wilson- Rae et al. (2007); Marquardt et al. (2007)) 
and technical ( [Schliesser et al.| ( |2008bD ; |Di6si| ( |2008l ); |Rabl et al.| ( |2009D ) lim- 
itations of radiation-pressure cooling of mechanical oscillators. In essence, if 
all technical sources of heating are avoided, the quantum fluctuations of the 
cooling light fleld provide a fluctuating force, driving the oscillator to random 
motion and therefore compete with the laser cooling eflFect. An advantageous 
ratio of these two eflFects can however be achieved in the so-called resolved- 
sideband (RSB) regime, in which the mechanical oscillation frequency ex- 
ceeds the cavity linewidth hi. In this case, cooling of the mechanical oscillator 
to its quantum mechanical ground state is possible. In the following, a brief 
outline of the required theoretical considerations will be given after a moti- 
vation by well-known results from atomic physics. In the main part of this 
section, experimental results — including the flrst demonstration of resolved 
sideband laser cooling of an optomechanical device — will be presented. 



6. 1 Ground state cooling: the atomic physics case 



The quantum mechanical expectation value of the energy of a harmonic me- 
chanical oscillator of Vt^ and mass meff is given by 

(i^mech) = + \m,s^l Ix') = (156) 



= hn^ ((n) + ^) (157) 

where the phonon number operator h = b'^b is given by the creation and 
annihilation operators 

b^ = — *— (x- i^^] (158) 
with the so-caUed zero-point fluctuations 



which are of the order of 100 am for typical silica microtoroidal resonators 
(meff = 10 ng, r2m/27r = 40 MHz). 
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|0,n- 1 



Fig. 34. Optical sideband cooling. Laser photons of energy hjjj\ — h{(jJo — Q^n) 
preferentially induce "red-sideband" transitions |0,n) \l^n — 1) (blue arrow) 
if they are detuned from the carrier |0,n) (gray arrow) and blue sideband 

|0,n) |l,n+ 1) (red arrow) transitions. As a consequence, the phonon occupa- 
tion n of the mechanical oscillator is reduced when the photon gets absorbed by 
the atom, or, in the optomechanical case, by the cavity. Re-emission of the photon 
(wavy lines), on average, does not change the phonon occupation (neglecting recoil). 



so that a detuned laser provides cooling. Figure from ref. (Schliesser et al. (2008b)). 



The question raised by the cooling results presented in the previous section is 
whether it is possible to reduce (i^mech) to levels comparable with the ground- 
state energy In other words, is it possible to reach the quantum ground 
state, in which the occupation number (the number of excitation quanta, 
phonons) (n) reaches zero? In this case, deviations from the classical cool- 
ing behavior described in the previous section are clearly expected in order to 
prevent cooling to reach arbitrarily low energy states. 



To answer this question, it is instructive to consider the results obtained in 
the context of laser cooling of ions (or atoms) (Hansch and Schawlow ( 1975[ ); 
Wineland and Dehmelt| ( |1975[ ); |wrneland and Itano| ( |1979D ; |Stenholm| ( |1986D ). 
Trapped in a harmonic potential, these elementary particles constitute me- 
chanical oscillators as well, their eigenfrequency being given by the tight- 
ness of the trap. Lasers can be used to drive electronic transitions of energy 
hujQ and lifetime /^~^, and optical and mechanical degrees of freedom are cou- 
pled by the Doppler shift of the optical resonance if the ion is moving and by 
the momentum transfer of absorbed and emitted photons. 



If the laser is detuned from the optical resonance, cooling can occur by favoring 
the absorption of a photon only in conjunction with the annihilation of a 
mechanical excitation quantum. Essentially, it is the phonon energy which 
makes up for the energy mismatch h{(jjQ — of the incoming photon to drive 
the electronic transition. If absorption takes place, the subsequently re-emitted 
photon has an average energy of hujQ (neglecting recoil). It therefore carries 
away the additional energy of the phonon, and leaves the mechanical oscillator 



in a state of lower excitation n (figure 34) 
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This method is usuaUy referred to as optical sideband coohng: an ion osciUat- 
ing in its trap exhibits absorption sidebands at frequencies ujq ± jQ^^ j G N, 



very similar to the case of a cavity discussed in section [271} If the laser is tuned 
to the red sideband at cjq — f^m, cooling transitions are resonantly enhanced 
as shown in figure 34 The quantum theory of laser cooling of trapped atoms 
or ions reveals that this method allows ground state cooling (n) pro- 
vided that ^ ^ (neglecting recoil) ( |Wineland and Itano ( 1979D ; Stenholm 



( [19861 )). In this case, the lowest average occupation that can be achieved is 
given by 

M = T7^- (160) 



Interestingly, this result can also be viewed as being due to the competition 
of the cooling eflFect of the detuned laser, due to an effective viscous force, and 
a "heating" effect due to quantum fiuctuations of the light beam, giving rise 
to a fiuctuating radiation-pressure force dltano et al.| (|1992D). 



6.2 Limits of radiation-pressure cooling using dynamical backaction 



A similar analysis can be applied to cooling by dynamical backaction. To assess 
the fundamental cooling limits, we consider the radiation pressure force as it 
is obtained from the quantum Langevin equations: 



-i{A + n) + k/2 +i{A-n) + K/2 



-i{A + n) + k/2 +i(A - + k/2^ 



(161) 



In essence, the expression in the first line correspond to quantum backaction 
due to quantum fiuctuations of the intracavity photon number, while the ex- 
pression on the second line is due to dynamical backaction (proportional to the 
mechanical displacement). The latter force contribution, is usually absorbed 



into the effective susceptibility Xefr from equation (27), while the spectrum of 
the force quantum fiuctuations is calculated as 



a;|pF2 V(A + J1)2 + (k/2)2 (A - 17)2 + (k/2)^ 



(163) 



72 



using the known correlation functions for the input quantum noise. We have 
introduced here the rates A_ and which can be shown to correspond 
to rates of ant i- Stokes and Stokes scattering events in which phonons are 



annihilated or created ( Wilson- Rae et al. (2007); Marquardt et al. (2007); 



Genes et al. (2008); Clerk et al. (2008)). 



The remaining energy stored in the mode after cavity cooling can be evaluated 
by integrating the mechanical displacement spectrum. This is possible analyt- 
ically if a hot mechanical reservoir hQ^ <C /cbT is assumed and the effective 
susceptibility still corresponds to a high-quality Lorentzian. Without further 
proof we note here that this can be safely assumed as long as the magnitude 
of the optomechanical coupling parameter 

G = 2a^o^zPF (164) 
is much less than For higher values of G, hybridization of optical and 



mechanical modes set in, a regime treated in detail in (Dobrindt et al. (2008); 



Groblacher et al. (2009a)). In the limits tv ^ G^T^ relevant to the work 
presented in the following, one eventually obtains (Genes et al. ( 2008] )) 



^m,snl{5x') = I Ixe^m' {S%i^) + (0)) 



dfl 
2^ 



1 



Teff 

r 



T^kBT + 



A. 

Teff ^ Teff 2 

For significant cooHng, one has A_ ^ A^ ^ Fm, and finally 



in) 




(165) 
(166) 
(167) 

(168) 



As a consequence, even for very strong cooling Feff ^ Fm, the phonon oc- 
cupation cannot be reduced to arbitrarily low numbers; instead it is bound 
by rimin- Two simple limits are derived depending on the ratio of mechanical 
resonance frequency and optical cavity linewidth: 

^min ^ jt;— ^ ^ uurcsolvcd sidebands (169) 

^min ^ — t;^:^ ^ rcsolvcd sidebands. (1^0) 

^ m 

Evidently, the ground state (n) can only be reached in the resolved- 
sideband regime, where <C f^m- Due to the non-zero occupation nbath ^ 
k^T /hTt^ when in equilibrium with the thermal bath, more stringent require- 
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ments arise in laboratory experiments. In particular, both k > Fmnbath (Do- 
brindt et al. (2008)) and > Fmnbath (Marquardt et al. ( 2007] )) are necessary 



to enable sufficiently high cooling rates without being limited by the cavity 
decay rate or an effectively overdamped mechanical resonator, respectively. 
Apart from these limits, it is desirable to work with optical resonances as 
narrow as possible to enable ground state cooling. 



Another advantage of working in the resolved-sideband regime is that the 
cooling rate increases monotonically for decreasing if all other parameters are 
fixed. As shown in figure [35| if the detuning is chosen for an optimum cooling 
rate, the intracavity power reduces as the laser is effectively detuned further 
from the WGM cavity resonance. This is advantageous to keep undesired 
nonlinearities at a low level, but also to prevent heating effects due to absorbed 
light. This effect is particularly relevant in cryogenic environments, where 
small absorbed powers may already significantly alter the temperature of the 



system, see section 7.3 



Finally, another technical imperfection in this cooling scheme may lead to 
a limit in the attainable occupation number. If the cooling laser frequency 
exhibits technical frequency noise, it is translated by the cavity into radiation 
pressure force noise. Laser frequency fiuctuations Sujoui^)^ with S^^{^) = 1 + 
4|5inp5'^a;(^)/^^ lead to additional radiation-pressure force fiuctuations with 
the spectrum 



^ a go 



(A - Q)2 + («;/2)2 {A + n)^ + {n/2)\ 



(171) 

which can be derived from the intracavity fiuctuations of the p-quadrature. 
A lower limit for the occupation can then be estimated by integration of 
the resulting displacement noise spectrum, taking into account the modified 
susceptibility of the mechanical oscillator (in particular its damping F^ba ^ 
Feff, cf. equation (29). In the linear cooling regime dPobrindt eFaL| ( |2008l )), 
the integrand is simply a Lorentzian, and we assume that the force noise is 
approximately constant within the frequency band in which it interacts with 
the mechanical oscillator. Then one obtains 



^FF (^m) 
2meffFdba 



(172) 



rifn 



As will be discussed later, this result gets modified by the additional presence 
of thermal noise (Schliesser et al. (2008b])), not taken into account in other 
discussions of frequency-noise induced cooling limits which (Diosi (2008); Rabl 
eFal] ( [20091 )). 
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Fig. 35. (a) Normalized cooling rate as a function of detuning and linewidth, for a 
fixed mechanical resonance frequency fl^ assuming G < k. Blue line indicates the 
optimum detuning for a given cavity linewidth. (b) Normalized cooling rate (blue 
line) and intracavity power (red line) as a function of linewidth for fixed if the 
detuning is kept at the optimum value. 

6.3 Resolved-sideband cooling of a silica microtoroid 



Implementing an optomechanical system which operates in the resolved-sideband 
regime is advantageous to suppress the effects of nonlinearities and absorption- 
induced heating present in most laboratory experiments. Furthermore, the 
analysis at the quantum level has shown that this regime is mandatory for 
cooling to the quantum ground state. It is interesting to note that this finding 
is quite generally valid, and applies in an analogous manner to a large class 
of systems in which a mechanical oscillator is cooled by parametric coupling 
to physical systems as diverse as an ion (Tian and ZoUer ( 2004D), a quantum 
dot ( Wilson- Rae et al. ( 2004D ), a Josephson qubit ( Martin et al.| (2004])), a 
superconducting quantum interference device (SQUID) ([Blencowe and Buks 



(2007)), a superconducting single-electron transistor (sSET) (Blencowe et al. 
j2005t ; |Naik et al.| ( [20061 )), ^ resonant electric circuit (Bro wn et al.| ( |2007D ). 



he excitations of all these systems need to have a lifetime exceeding the 
oscillation period ^ of the mechanical oscillator to enable ground-state 
cooling. 



Experimental implementation of the resolved-sideband regime however proves 
non-trivial: in a classical optomechanical system involving a Fabry-Perot res- 
onator, increasing the cavity length L in principle reduces but at the same 
time also the coupling strength go = —uj^^/L. Also, it typically increases the 
waist size of the beam, which may induce diffraction losses (and therefore in- 
crease K again) if, as it is typically the case, the micromechanical oscillator 
constitutes one of the cavity mirrors. In an integrated system such as silica 
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microtoroids, increasing the cavity size (within the relevant range of 60 to 100 
/im diameter) typicaUy leads only to modest improvements of as the losses 
are dominated by absorption along the light propagation path. In addition to 
the reduced coupling go = —uj^^/R^ the mechanical resonance frequency of the 
RBM is reduced for larger cavity size. The first experiment demonstrating ef- 
ficient resolved-sideband laser cooling was eventually performed in 2008 with 



silica microresonators (Schliesser et al. (2008b)), while other systems, such 



as superconducting microwave cavities ( [Teufel et al. (2008)) or Fabry- Perot 
resonators (Groblacher et al. (2009b)) followed soon thereafter. 



Figure [36] shows data obtained with a 47 /xm-diameter toroid hosting a RBM at 
^m/27r = 73.5 MHz. A first laser drives the RBM to regenerative oscillations, 
by pumping a WGM in a blue-detuned manner, leading to the well-known os- 



cillatory instability 


(Kippenberg et al. ( 


2005D;|Carmon et al. 


(2005 


); 


Rokhsari 


et al. 


(2005, 


2006 


);] 


iossein-Zadeh et al. 


(2006|)). A second laser, launched into 



the same tapered fibre, is used to probe a diflFerent, high-Q WGM resonance 
of the oscillating toroid. Using optical filters to extract only the transmission 



signal of the probe laser we obtained the traces shown in figure [36] as this laser 
is swept over the resonance. 



13=0.94 




P=1.75 



-200 -100 +100 
Detuning (MHz) 



+200 



Fig. 36. Steady- state transmission spectrum of a microtoroid when the mechanical 
degree of freedom is excited to picometer-scale amplitude oscillations at its reso- 
nance frequency of Vt^/2Ti — 73.5 MHz using an auxiliary laser. The linewidth of 
the resonances corresponds to the optical decay rate of A^/27r = 3.2 MHz, placing 
this device deeply into the resolved sideband regime. Blue points are experimental 



data, and the red lines are fits using the Bessel function expansion (equation (173 



with /3-parameters given in each panel. Figure from ref. (Schliesser et al.|(2008b| 



Instead of a single dip, several dips are observed if the electronic signal is 
filtered with a low-pass filter with a cutoff far below f^m- From the intracav- 
ity mode amplitude of the oscillating cavity (equation (jTl)) the low-frequency 



76 



Monitoring 

l^ser Monitoring 
cavity 
mode 



Cooling 
laser 





Cooling 
/ cavity 
mode 



Optical 
frequency 
► 



Fig. 37. Scheme used for resolved-sideband cooling. The cooling laser at A ^ 980 nm 
(green line) is tuned to the lower sideband (A = — f^m) of a high-Q optical resonance. 
Resonant anti-Stokes scattering into the cavity mode, at rate A_ dominates over 
Stokes scattering at rate A^. A second laser is tuned in resonance with a different 
WGM at A = 1064 nm, and used to monitor the mechanical displacements. Figure 



from ref. (Schliesser et al. (2008b)). 



component of the transmitted power 



-DC 1 2 
^out I 



1 -^c(l -^c)^ 



(173) 



is calculated ( [Schliesser et al. (2008b)), where are the Bessel functions 
and /3 = S'o^o/^m- From fits using this model, it is possible to derive the 
amplitude Xq of the mechanical oscillation as about 5.3, 8.3 and 9.9 pm for 
the three traces, respectively. At the same time, these traces clearly show that 
the mechanical resonance frequency — by which the Lorentzian dips are 
spaced — largely exceeds the optical cavity linewidth tv. Indeed, the fits yield 
a resolved-sideband factor Q^/tv ^ 23, due to the very high cavity finesse of 
440,000. 

For resolved-sideband cooling, the two lasers are used in a different way. 
The cooling laser is tuned below the narrow WGM, to the lower sideband at 
A = — f^m- The second laser is used to sensitively monitor mechanical displace- 



ments and therefore is tuned in resonance with a different WGM (figure 37). 



Using two completely independent laser sources helps to rule out artifacts in 
the displacement measurements, due to potentially induced noise correlations 
between the cooling beam and induced motion of the mechanical oscillator, 
referred to as "squashing" in feedback cooling (Poggio et al. ( 2007] )). 



In this experiment, sensitive monitoring of mechanical motion is accomplished 
using the Hansch-Couillaud technique described in section |4.2.2[ As shown in 
figure |38| a low-noise Nd:YAG laser is used for this purpose, which is locked to 
the line-center using feedback to a piezo in contact with the YAG crystal. In 
order to lock the cooling laser to a detuning much larger than the resonance 
linewidth we use an experimental technique, in which the signal obtained 



by frequency modulation spectroscopy (Bjorklund et al. (1983)) is used as an 
error signal. 
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Fig. 38. Setup used to implement resolved-sideband cooling and high-sensitivity 
monitoring of the RBM of silica toroidal microresonators. Two lasers are used, a 
diode laser at 980 nm for cooling by tuning and locking it to the lower sideband of a 
high-Q WGM, and a Nd: YAG laser at 1064 nm monitoring mechanical displacements 
using the Hansch-Couillaud technique. More details are given in the text. Figure 



from ref. (Schliesser et al. (2008b)) 
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Fig. 39. Displacement noise spectra of the RBM as recorded by the monitoring 
laser with the cooling laser off (left panel), and running at a power of 2.7 mW, 
when detuned to the lower sideband (right panel). Cooling rates up to 1.6MHz 
can be achieved. Circles are data points and lines Lorentzian fits. Figure from ref. 
dSchliesser et al.| ( |2008b"l )). 



Figure [39] shows cooling results obtained with the 73.5 MHz-sample. Note that 
due to the strongly detuned operation, only a fraction of about ^ (4(f]ni/^)^ + 
1)~^ 5 X 10~^ of the launched power of 3mW, i. e. circa 1.5 /iW, is coupled 
into the cavity. Nonetheless, very high cooling rates up to Fdba/STT = 1.56 MHz 
can be achieved. The mode temperatures achieved with this sample, however, 
remained above 20 K, due to the modest mechanical quality factor of about 
= 2,100, and a heating mechanism described below. 
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Fig. 40. Resonance doublet observed in the power spectrum of the coohng laser for 
strong cooling powers. See text for more information. 

We have also recorded noise spectra of the mechanical oscillator by directly 
analyzing the transmitted power of the cooling laser (as explained in section 
[5]). Qualitatively the same behavior of the mechanical spectra is observed 
during cooling, however with significantly worse signal-to-noise ratio. This 
renders, for example, the determination of the mode temperature extremely 
difficult for strong cooling. Furthermore, for cooling laser powers on the order 
of 3mW, significant deformations of the mechanical spectra from their orig- 
inally Lorentzian shape were observed. The particular shapes were strongly 
dependent on the detuning and coupling conditions. Figure |40l shows one ex- 
ample in which a resonance doublet appeared. We note that such an doublet 
would in principle be expected when entering the regime of strong optome- 



chanical coupling, in which optical and mechanical modes hybridize (Dobrindt 
e Fal] ( [2008| ; |Groblacher"eraL] ( [2009a| )). 



To achieve lower mode temperatures, a second cooling run was initiated. We 
used a larger toroid (radius R = 28 /xm) with a RBM at Q^/27t = 40.6 MHz 
and r^/27T = 1.3 kHz, corresponding to a very high mechanical quality factor 
of 30,000. The broadband displacement spectrum recorded with this sample 
the Hansch-Couillaud technique (cooling laser is oS) is shown in fig- 



usme 



ure kll A displacement sensitivity at the 10~^^ m/yHz-level is achieved. At 
low Fourier frequencies, the noise spectrum again reveals a thermorefr active 
background already discussed in section |4} 

The strong peaks observed at 14.96, 28.58 and 40.59 MHz can be assigned to 
different radially symmetric modes in this sample. The latter peak belongs to 
the RBM, which has an eflFective mass of about 10 ng. If the cooling laser is 
tuned to the lower sideband of a /^/27r = 5.8 MHz- wide WGM, cooling of the 



RBM is evident by the reduction of the RBM's thermal noise (figure 41). It 



is interesting to note that the thermal noise in the other radially symmetric 
modes is not aflFected. In the resolved sideband regime, this is possible due to 
the fact that the absorption sidebands of the individual mechanical modes do 
not overlap, as they are as narrow as the optical resonance. Furthermore, we 
note that eflFects analogous to "sympathetic cooling" — due to, for example a 
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Fig. 41. Resolved-sideband cooling of the RBM of a microtoroidal oscillator. 
Top panel shows a broadband displacement noise spectrum recorded using the 
Hansch-Couillaud technique (red trace). Gray trace is background trace with the 
fiber taper retracted from the WGM near field. It is dominated by shot noise, its 
dependence on Fourier frequency is due to the weaker transduction of displacement 
fluctuations to phase fluctuations at Fourier frequencies beyond the cavity cutoff 
The sensitivity is at the lO"-*^^ m/\/Hz-level. Insets show the displacement pat- 
terns of three radially symmetric modes, with the corresponding strain indicated in 
the color code. The bottom panel shows zooms on the displacement noise of these 
three modes, when the cooling laser is off (red traces) and on (blue traces), at a 
detuning of A = — f^m- Clearly, only the RBM mode (number 3) is affected. The 
lowest achieved occupation number in these measurements was (n) ^ 5900. Figure 



from ref. (Schliesser et al. (2008b)) 



preferential coupling of the RBM to other radially symmetric modes- 
observed. 



-is not 



The mode temperature of the RBM is subsequently evaluated by integrat- 
ing the displacement noise spectrum Sxx{^) (cf. section 5.2). The lowest mode 
temperature obtained with this sample is ^ 11 K, corresponding to a resid- 
ual occupation of (n) ^ 5900. As in the experiments with the first sample, this 
is significantly above the value expected from the high cooling rates, reaching 
up to Fdba/STT = 119 kHz with this sample. This discrepancy is attributed to 
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excess heating by the coohng laser's frequency noise. 



Indeed, an independent measurement reveals frequency fluctuations of the 
cooling laser at the level of 



^Sujuj{Qrn) = 27r • 200Hz/vHz, a value consistent 
with earlier measurements on the frequency noise at radio frequencies of a 



grating stabilized diode laser (Zhang et al. ( 1995| )). For a flnite bath temper- 
ature T, the lowest occupation that can be attained can be estimated by 



(n) 



2meffrdba 



(174) 



In the resolved-sideband limit <C f^m, the force noise due to frequency 



fluctuations from equation (171) simplifles to 



^FF (^m) 



(175) 



and the cooling rate is approximately given by 



dba 



K 



(176) 



As a consequence, the lowest temperature is obtain ed at a flnite input powe r, 
and one flnds nmin ' 



^2kBTmeET^S^^{n,n)/h\go\ (Schliesser et al. (2008b)). 
For the parameters of this experimental run, T = 300 K, mgff = 10 ng, i m/27r = 

1.3kHz, ^Js^u;{^m)/^l, ^ 4/xrad/yife, R = 38 /xm, n^/27r = 40.6 MHz, 
u)/27r ^ 300 THz, one obtains ^ 5200, close to the lowest observed occu- 
pation number. 



6.4 Direct sideband spectroscopy 



If the ground state of the mechanical oscillator is approached, it becomes in- 
creasingly diflicult to reliably measure its occupation number. For example, 
if the method demonstrated here — using an independent motion transducer 
based on an optical phase measurement — the requirements on its performance 
are quite challenging. Even for a perfect implementation of such a measure- 
ment, in which the standard quantum limit can be reached (cf. section [4]), the 
signal-to-background ratio at the mechanical resonance frequency would be 
only equal to unity at the mechanical resonance frequency when the resonator 
is in its ground state: for (n) 0, the spectrum of displacement fluctuations 
is expected to peak at a value of 
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which exactly equals S^^^{Qm)' 

We note that our experiment reaches an imprecision level close to the SQL: For 



example, for the sample described above, we calculate J S^^{Qrn) ^ 2.2 am/v Hz 



and routinely achieve imprecision noise at the same level (cf. figure 41). How- 
ever, as in most other experiments reporting similarly low imprecision levels 
(Regal et al. (2008); Teufel et al. ( 2009] )), backaction of these measurements 



is diflScult to quantify and may impede cooling to the quantum ground state. 

As an alternative approach to measure the residual occupation number for 
small (n), it has been suggested to directly monitor the individual motional 
sidebands of the cooling laser ( Wilson- Rae et al. ( 2007[ )). In theory, the optical 



spectrum of the cooling laser light coupled back into the taper is described 
by the expression 



+ 



rye A2 + (/,/2)2 

reff/(2^) 




(^_(^^ + ^J)2 + (reff/2)2 

reff/(27r) 

■(^_(^^_f]J)2 + (reff/2)2 




(178) 



where $ is the photon fiux emerging from the cavity ( [Wilson- Rae et al." ( 2007[ )). 



Next to the central line at the carrier frequency cji, two sidebands are ex- 
pected, their lineshape being determined by the eflFective susceptibility of the 
mechanical oscillator. Most notably, the spectrally integrated photon fiux of 
the individual sidebands is given by A_ (n) and ((n) + 1) for the upper 
and lower sidebands, respectively. As a consequence, the sideband asymme- 
try, given initially by the asymmetry in A_ and ^4+, becomes balanced for 
suflBciently low (n) . The change in the ratio of the sideband amplitudes could 
therefore serve as a gauge of the occupation (Wilson- Rae et al. ( 2007] )). 



Accessing the individual sidebands in a measurement necessitates the ability to 
individually resolve them against the much stronger carrier signal at the laser 
frequency. For typical parameters of an optomechanical experiment, this is dif- 
ficult, as in the RSB regime the power even in the resonantly enhanced upper 
sideband is weaker by ^ (n)5'o4r7^x|pp/f]^ ^ {n)O{10~^) than the carrier, and 
the lower sideband is again weaker by a factor of /^^/16f]^. As a consequence, 
sufiicient suppression of the carrier using a (single) filtering cavity is hardly 
possible, as it would require sub-kHz cavity linewidth, which could not even 



Note that here the high-Q approximation for the mechanical susceptibility was 
used, ((J^^-f^2)2_j^2p2^)-i ^ (4^^2^((f7-J^^)2 + (reff/2)2))-i. Wealso emphasize 
that this spectrum lies in the optical domain, instead of the RF-domain of the 
spectra discussed previously. 
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simultaneously collect all light in the sidebands (typically reff/27r > IkHz). 



An alternative way of individually resolving the sidebands is to use a hetero- 
dyne technique. A similar technique has been demonstrated to enable resolving 



motional sidebands of a laser-cooled ion (Raab et al. ( 2000| )). In such an 
experiment, the cooling light which couples back to the fiber taper is mixed 
with a strong local oscillator beam at a diflFerent frequency uj\ — Vt^oM with 
^AOM > ^m- Then the upper sideband, carrier, and lower sideband signals 
are detected at the radio frequencies J^aom + ^m, ^aom and J^aom — ^m, 
respectively. 



Figure |42] shows the implementation and results of a proof-of-principle exper- 
iment using this technique. When the laser is tuned close to resonance, the 
motional sidebands have roughly the same amplitude, which are again inde- 
pendently calibrated using a frequency-modulation technique. As (n) 1, 
the residual asymmetry in this case is attributed to a finite detuning from the 
WGM resonance. When the laser is further detuned, the asymmetry becomes 
more pronounced, as A_ ^ In this experiment, an asymmetry of more 
than 15 dB is reached. Higher values could not be observed due to the limited 
SBR. This limitation is due to the classical frequency noise of the cooling laser, 
which could be shown to induce the relatively high background noise level in 
this measurement. 

We note however that, in contrast to the experiment performed on trapped 
ions (Raab et al. ( 2000[ )), this experiment was performed without actively 



driving the motion of the mechanical oscillator. In our case, the oscillator is 
exclusively driven by thermal noise. In that sense, this measurement tech- 
nique beautifully reveals in a very direct manner the physical process under- 
lying cooling by dynamical backaction: enhanced anti-Stokes scattering and 
suppressed Stokes-scattering. 



Note however that only micro-motion sidebands could be detected for a cold ion. 
To observe the secular motion, the ion was driven to large-amplitude oscillations in 
the trap. 
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Fig. 42. Heterodyne spectroscopy of the motional sidebands during cooling, (a) 
Experimental setup, in which a frequency-shifted laser beam is used as a heterodyne 
local oscillator, (b) Experimental results, showing the power spectral density of 
the differential photocurrent recorded in the heterodyne receiver. Red points are 
recorded with the cooling laser close to resonance, and blue points are recorded 
with a detuned coohng laser. Figure from ref. (Schliesser et al. ( 2008b[ )). 
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7 Approaching the quantum ground state 



The work presented in the previous sections has clearly demonstrated the po- 
tential of resolved-sideband cooling. In the following, we present experimental 
results in which this technique is successfully implemented to cool a mechani- 
cal oscillator close to its quantum ground state, rivalling the results achieved 
with nanoelectromechanical systems engineered for this purpose for nearly a 



decade (Knobel and Cleland (2003); LaHaye et al. (2004); Schwab and Roukes 



(2005); Naik et al. (2006); ?). In our laboratory these results (Schliesser et al. 



(2009)) were achieved by combining laser cooling with efficient cryogenic pre- 
cooling, a technique traditionally applied to nanoelectromechanical systems, 
but novel in the context of optomechanics. 



As a second crucial feature, the vastly superior displacement sensitivity of 
our optomechanical system warrants the ability to monitor such an ultracold 
oscillator at levels close to the amplitude of its zero-point fluctuations Xzff = 
^Jh/2mQf[Q^ — in spite of the fact that its mass is more than 1000-times higher 
than in typical nanomechanical systems. 

Finally, considered from the perspective of quantum measurement, we are 
able to assess the backaction of the measurement performed on the ultracold 
oscillator. We will show that considering both imprecision and backaction of 
our measurement, one can conclude that the optical displacement transduction 
performs in a near-ideal manner. 



7. 1 Implementation 



From the theoretical considerations of the previous sections, we can summarize 
the limits in cooling by dynamical backaction in the simple expression 

Reaching (n) therefore requires low mechanical dissipation Fm, high cool- 
ing rates F^ba? a low environment temperature T, high mechanical resonance 
frequency low frequency noise S^^{Qrn) ^ ^fn of the cooling laser (or, 
more general, no technical noise of any kind), and operation in the resolved- 
sideband regime to achieve low A^/{A_ — A^) ^ /^^/16f^^. 

After the optimization of the silica microtoroids for low dissipation, high fre- 
quency and operation deeply in the resolved-sideband regime, and the demon- 
stration of very high cooling rates, the most signiflcant performance improve- 
ments are expected for operation in a cryogenic experiment, and using a cool- 
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ing laser devoid of any technical noise at the relevant radio frequencies. In the 
following, we briefly describe the implementation of these two advances. 



7.1.1 Cryogenic optomechanical experiment 



For the successful cryogenic operation of an optomechanical cooling exper- 
iment, it is important that the concrete technical implementation provides 
a proper thermalization of the sample to the cryogen, even in presence of 
a microwatt-scale heat load through residual light absorption in the sample. 
This may be a non-trivial task, as strong thermal anchoring of the mechan- 
ically compliant device would in many cases imply designs opposite to the 
requirement of low clamping losses of the mechanical oscillator. For example, 
the highest mechanical quality factors in silica microresonators were achieved 
with toroids held only on a sub-micrometric pillar (Schliesser et al. ( 2008b[ )), 
or suspended from silica spokes (Anetsberger et al. (2008)). Heat transport 
from the region of the optical mode in the torus to a heat sink below the 
silicon chip can be expected to be prohibitively low. 



For this reason, we chose to directly immerse the sample into the cryogen, 
a dilute helium gas, similar to an early experiment with silica microspheres 
(Treussart et al. ( 1998[ )) immersed in superfluid helium (flgure 43). In con- 
trast to other cryogenic experiments with optical microcavities employing a 
cold flnger cryostat ( ^Srinivasan and Painter (2007); Groblacher et al. ([2008')), 
this allows us to reliably thermalize the sample to the base temperature of 
the cryostat (1.6 K). At the same time, due its low pressure (typically held 
in the range 0.1 mbar-50mbar), and the increased intrinsic damping of the 
mechanical mode due to two-level fluctuators (see section 3.2.4), the helium 
gas present in the sample chamber does not induce signiflcant damping of the 
mechanical mode. 



The employed commercial Optistat SXM cryostat (Oxford instruments) im- 
plements such a cooling scheme by providing two completely independent 



helium containers (cf. flgure 43): A larger (volume 4.3 liters), thermally in- 
sulated reservoir contains liquid helium. From this reservoir, a rotary pump 
continuously draws helium through a capillary contained in a copper block. 
During this process, the helium evaporates and cools to temperatures down 
to 1.6 K. The copper block has a large central bore (49mm diameter), and 
constitutes part of the wall of the second helium container — an approximately 
50 cm-long tube flUed with low-pressure helium gas. It therefore serves as a 
heat exchanger between the cold He gas from the reservoir and the buffer 
gas in the central chamber in which the sample is held. The temperature of 
the heat exchanger can be continuously tuned by controlling the evaporation 
rate of the He using a control valve. In addition, an electric heater in the 
copper block allows electronic temperature stabilization, and makes the whole 
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Fig. 43. Implementation of cryogenic cooling. Left panel: schematic drawing of the 
employed buffer gas cryostat. Liquid ^He, kept in an isolated reservoir, evaporates 
when pumped through a capillary ("heat exchanger") and thereby cools low-pres- 
sure (0.1 mbar-50 mbar) ^He-gas in the sample space to temperatures down to 1.6 K. 
However, a heater can also be used to stabilize the temperature of the buffer gas to a 
different value in the full range up to 300 K. The experimental assembly ("head") is 
directly immersed into the buffer gas. Right panel: side (top) and bottom (bottom) 
view of the experimental head. The coupling taper, glued to a glass taper mount, 
is fixed, while the position of the chip containing the silica microresonators can be 
adjusted using stick-slip piezoelectric actuators. Optical access from two orthog- 
onal directions allows convenient addressing the toroids on a chip, and precision 
positioning for coupling to the tapered fiber. 



temperature range up to 300 K accessible. 



Successful coupling to the WGM in the toroids requires a stable, micrometer- 
scale gap between the fiber taper and the edge of the toroid. Low vibra- 
tions in the experimental assembly are therefore an important selection cri- 
terion. The chosen system has successfully been operated in many labs for 
the purpose of low-temperature scanning-probe microscopy — obviously very 
vibration-sensitive applications. Both the suspension of the experimental head 
from a ^ 0.5m-long metallic rod into the buffer gas, and our compact head 
design (figure 43) indeed renders the cryogenic coupling setup extremely stable 
(Riviere et al. (2010)). Standard techniques to reduce the coupling of vibra- 
tions from the running pumps (heavy masses on the transfer pipes between 
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cryostat and pump) are sufficient to eliminate vibrations to a degree that they 
cannot be observed in any experimental signal. At the same time, this cryo- 
stat allows convenient optical access from two orthogonal directions, providing 
both a side and top view of the silicon chip. Standard machine vision prod- 
ucts (Navitar 12 x) with large working distance 9 cm) allow observation 
of the coupling region with up to x 14-magnification, sufficient for convenient 
addressing and alignment of the toroids versus the fixed fiber taper. 



Preparation of a cryogenic cooling experiment includes two critical tasks: se- 
lection of an adequate cooling sample and preparation of a fiber taper mounted 
into the experiment. In particular, the resonators on the chip are tested for 
high-Q WGM resonances allowing to enter the RSB regime (/^ < Q^) and for 
low-mass, high-quality RBM modes. These measurements are done in stan- 
dard room-temperature coupling setups. If an adequate sample is found, it is 
mounted into the cryogenic head. 



The compactness of the cryogenic experiment renders the use of the bulky 
metallic taper holders used during taper fabrication impossible. Instead, tapers 
have to be transferred to a compact glass support made of a microscopy glass 
slide which, importantly, features the same coefficient of thermal expansion as 
the taper. In that manner, the taper tension remains constant during cooldown 
of the cryostat. 



The fabrication and transfer of the fiber taper proceeds in several steps which 
are illustrated in figure 44, After the standard fabrication procedure in a hy- 
drogen fiame, the taper, still in the fabrication holder, is placed in an auxiliary 
testing setup. Here, a microtoroid is approached until it touches the taper at 
its central position. If retracted again, the toroid pulls the taper with it for 
a certain distance due to proximity forces. Using a micrometer drive, the ta- 
per is strained by increasing the separation of the clamps ("elongation" in 
figure 44), and the tests are repeated until the distance over which the taper 
can be displaced by the toroid is in the range of 10 /xm. While not particularly 
quantitative, this simple test was found very helpful in the delicate task of 
mounting a well-strained taper into the cryostat: fioppy tapers make coupling 
at low temperatures impossible, while too tightly strained tapers are prone 
to rupture during cooldown. The glass support, fabricated from a simple mi- 
croscopy slide is then positioned underneath the taper with a 3d-translation 
stage. The taper is glued to the support with UV-hardening epoxy applied 
first on one side, hardened, and then on the other side. Then the fiber is re- 
leased from the metal holder and transferred to the cryostat. We finally note 
that due to its reduced diameter and potentially lower quality of the available 
single-mode fiber material, creating a cryogenic fiber taper at shorter wave- 
length (e. g. 780 nm as in this work) is significantly more difficult than at 
near-IR wavelengths (1.5 //m). 
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1. taper pulling (side view) 



2. taper tensioning (top view) 





3. positioning of cover slide (top view) 





4. taper glueing (top view) 




^ ^ 




5. removal of metal taper holder (top view) 
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Fig. 44. Five steps in the fabrication and transfer of a fiber taper to the compact 
glass support used in the cryogenic experiments (left panel). After pulling the taper 
over the flame of a hydrogen torch, its strain is tested using an auxiliary chip with 
microtoroids (upper right graph). If properly strained, the glass support (prepared 
from a cover slide) is positioned underneath the taper, and the fiber is glued to the 
support with an epoxy. After hardening the epoxy with UV light, the metal clamps 
of the fabrication holder are opened to release the taper now only held by the glass 
support (photograph bottom right). It can then be mounted into the experimental 
head. See text for more details. 

The slide carrying the taper is then mounted into the experimental head, and 
the fiber ends are guided to the top of the sample insert and leave the cryostat 
via a helium-tight feedthrough. We have achieved total optical transmission 
through the cryostat up to 80% from fiber end to fiber end. A toroid from 
the mounted sample can then be approached using the piezoelectric stick- 
slip actuators (Attocube systems) that carry the sample chip. For testing 
purposes, and to check the correct alignment of the taper mount (position, 
tilt), these experiments can also be done with the sample insert (experimental 
head, suspension rod, and top vacuum flange) held outside the cryostat. 

To verify the effectiveness of the buffer gas, a sample with microtoroids, and 
a coupling taper suited for 1550nm-light was mounted into the cryostat. The 
cryostat was cooled down, and the mode temperature of this sample's RBM 
was determined using a displacement measurement based on the PDH tech- 



nique (see section 4.2.3). To enhance the sensitivity, the weak probing light 
(< 2 /iW) from the employed low-noise fiber laser (BASIK, Koheras) was am- 
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plified using an erbium-doped amplifier after leaving the cryostat, allowing us 
to measure displacement spectra down to the base temperature of the cryostat 
(Arcizet et al. ( 2009[ )). In figure 45, the mode temperature extracted from 
the calibrated spectra is shown to closely follow the temperature of the helium 
buffer gas down to about 1.8 K, or (n) ^ 600. In spite of being exclusively in 
contact with the cold helium gas, this 62MHz-RBM could therefore be ther- 
malized to occupation numbers which are, for lower frequency oscillators, often 
only attainable in dilution refrigerator systems. 




1 10 

Cryostat temperature (K) 

Fig. 45. Thermalization of the 62MIIz-RBM of a silica microtoroid. The mode 
temperature Tm follows the temperature of the helium buffer gas down to 1.8 K, 
or (n) ^ 600, proving the effectiveness of the buffer gas cooling. Figure from ref. 
dSchliesser eFall ( [20091 )). 



Cryogenic cooling of the silica microresonators implies severe modifications of 



the properties of these devices. As discussed already in section [3. 2. 4[ the me- 
chanical quality factor critically depends on the temperature of the structure, 
due to coupling of the mechanical modes to two-level systems in the glass. In 
contrast, the change in the sound velocity, and therefore the mechanical reso- 
nance frequencies, is a relatively weak effect at the (9(10~^)-level ( [Arcizet et al. 
( [20091 )). At the same time, the thermorefractive (dn/dT) and thermoelastic 



(a) coefficients are strongly modified. While a tends to zero at low tempera- 
tures (T < 2 • 10~^ at 1.6 K), the thermorefractive coefficient continuously 
drops and changes sign at a temperature of about 8K. The resulting non- 
trivial temperature dependence of the optical resonance frequency gives rise to 
interesting thermal multistability effects when an optical resonance is probed 



using a laser powerful enough to heat the cavity above this temperature (Ar- 



cizet et al. ( 2009D ). In the context of the cooling experiments, it is important 
to note that negative dn/dT implies that the red wing of an optical resonance 
is dynamically self-stable, relaxing dramatically the requirements on the laser 
frequency stability and locking speed when working with red detunings. 



13 - 



(Treussart et al. (1998)). 



Note these observations (Arcizet et al. (2009)) differ from a previous report 
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7.1.2 Optical system 



For a cooling experiment with the aim of reaching very low occupation num- 
bers, it is mandatory to avoid excess backaction noise in the cooling laser 



beam, as this leads to a limit ([Schhesser et al. (2008b); Diosi (2008); Rabl 



et al. ( 2009| )) in the cooling performance. The grating-stabilized diode lasers 



employed in our previous cooling work was observed to exhibit excess fre- 
quency noise (cf. section 6.3). One possible strategy to eliminate this noise is 
to measure the noise using an independent, quiet cavity, and fast feedback on 
the frequency of the laser using e. g. the diode laser current as an actuator. 
Feedback with the required bandwidth exceeding 50 MHz however is techni- 
cally difficult, as a few meters of cable alone would already induce phase lags 
in the loop sufficient to render it ineffective. 



As an alternative approach, a quiet cavity could be used in transmission to 
filter out the noise. In this manner, reduction of phase noise by 85 dB has 
been demonstrated at IMHz Fourier frequency (Hald and Ruseva ( 2005] )). 
However, the filtering reduces the available power by more than a factor of 5, 
so that the remaining available power of ^ 2mW is problematically low for 
cooling experiments. At the same time, the complex filtering setup derogates 
one important asset of a grating-stabilized diode laser: its convenient, wide- 
range, and fast tunability. 



For these reasons, we chose to use a solid-state pumped continuous wave 
titanium-sapphire (Ti:S) laser (Matisse TX, Sirah Lasertechnik) , which com- 
bines very wide tunability (750 to 870 nm) with high output power (> 0.5 W) 
and quantum-limited noise in both amplitude and phase at the relevant Fourier 
frequencies. Figure [46] shows the schematic setup of this device. The laser is 
pumped with up to 6 W at 532 nm from a frequency-doubled diode-pumped 
solid-state laser (Millenia, Spectra Physics). The frequency of the laser is ad- 
justed by a combination of intracavity filters: a Lyot birefringent filter narrows 
down the frequency range in which gain exceeds loss to a few hundred GHz, 
and a thin and a thick etalon, with free-spectral ranges of about 250 GHz and 
20 GHz, respectively, single out one longitudinal mode of the laser to oscillate. 
Once the desired laser mode oscillates, electronic feedback loops are used to 
lock the etalons to keep their transmission maximum resonant with the laser 
mode. These loops are implemented by the controller of the laser and can be 
adjusted on a computer using an USB-interface. 



Two mirrors mounted on piezo-electric transducers can then be used to fine- 
adjust the frequency of the laser. In particular, the "slow" transducer allows 
frequency scans up to about 60 GHz. The "fast" transducer, together with the 
intracavity electro-optic modulator, can be used for fast feedback to the laser 
frequency. In normal operation, these actuators are employed in a feedback- 
loop to lock the laser to an external, temperature-stabilized, ^ 30 cm-long 



91 



slow 
piezoelectric 
transd^^^ 



pump beam 




thick etalon 



Faraday rotator 



fast 

piezoelectric 
transducer 



electo-optic 
crystal 



thin etalon 



y 

Lyot filter / output coupler 
out-of-plane mirror 



Fig. 46. Layout of the titanium-sapphire laser used for the cryogenic coohng exper- 
iments. The Ti:S crystal in the ring cavity is pumped optically at 532 nm (green 
beam). In order to select a single longitudinal mode within the tuning range of 
750 to 870 nm (determined by the cavity mirrors), a Lyot filter, thick etalon and 
thin etalon are used. A Faraday rotator, one reflection point lying outside the plane 
of the rest of the cavity mode, and a number of Brewster-angled surfaces ensure 
unidirectional lasing action. The precise frequency of the laser can be actuated by 
displacing cavity mirrors using a fast, a slow piezo, as well as an intracavity elec- 
tro-optic modulator. 

high-finesse {J-" ^ 300) cavity via the PDH technique, in order to eliminate 
low- frequency frequency noise of the laser, reducing the laser line width to 
< 30 kHz. In this case, the laser frequency is controlled by tuning the reso- 
nance frequency of the reference cavity, again by displacing a cavity mirror 
mounted on a piezoelectric transducer. The custom version of the laser set 
up in our laboratory allows alternatively to directly lock the laser to a WGM 
in a silica microtoroid, by feeding an error signal (derived using one of the 
above-described methods) to the laser controller. 



The Ti:S laser can be tuned to arbitrary frequencies within its tuning range; 
however, continuous scanning over the whole range is not possible, as the 
various frequency selectors (Lyot and etalon filters) have to be readjusted if 
the laser has to be tuned by more than about 30 GHz. We therefore resort to 
a continuously tunable grating-stabilized diode laser (DL) in Littman-Metcalf 
configuration (TLB-6300, New Focus) for pre-characterization of the toroids 



(cf. figure 47). Once a well-suited WGM resonance is found with the DL, the 



Ti:S laser is tuned close to the frequency of this resonance. 

For this purpose, the Ti:S frequency is coarsely adjusted first, using a grating 
spectrograph. For the subsequent fine-adjustment, it was found very useful to 
monitor the interference signal of the DL — scanning over a broader frequency 
range including the WGM resonance — and the slowly tuning Ti:S laser. When 
the DL scans through the Ti:S frequency, a transient interference signal is 
observed as long as the difference frequency is within the bandwidth of the 
receiver (125 MHz). This transient beat serves as a convenient marker of the 
Ti:S frequency, while the standard transmission signal of the DL indicates the 
WGM resonance frequency. 

After the Ti:S has been tuned close to the WGM resonance, an optical switch 
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Fig. 47. Schematic of the setup used for resolved-sideband laser cooling in a cryo- 
genic environment. An external cavity diode laser (ECDL) is used for precharac- 
terization of the optical modes. When a suited WGM resonance is found, the Ti:S 
laser is tuned to the frequency of the ECDL using the "diagnostics" optical out- 
put (see text). Then the input to the experiment is switched to the Ti:S using a 
MEMS switch. The subsequent optical setup is essentially a Mach-Zehnder inter- 
ferometer, one arm of which contains the cryostat with the taper coupled to the 
WGM in a microtoroid. For testing purposes, the output of the taper can also di- 
rectly be monitored by guiding the light to an auxiliary photoreceiver (dashed beam 
path). In normal operation, the transmitted light is spatially overlapped with the 
local oscillator in a polarizing beam splitter cube, albeit in orthogonal polarization 
modes. Projected into a different polarization basis using a waveplate and another 
polarizing cube, the modes interfere, and the interference signal is recorded with 
a balanced receiver. The phase of the local oscillator — and therefore also the de- 
tected signal quadrature — is locked using the DC-signal from the balanced receiver. 
The radio-frequency-part is sent to a spectrum analyzer to obtain the displacement 
noise spectra, calibrated again by the frequency modulation of probing laser using 
an elecro-optic modulator (EOM). Waveplates and polarizing beam splitter cubes 
are used at various positions in the optical setup to adjust powers and polarizations 
of the beams. 

based on MEMS -technology is used to switch to the Ti:S laser as input to 
the experiment. The optical setup of the experiment corresponds essentially 



to a homodyne measurement (figure 47). To achieve the highest possible SBR, 



good interference contrast between the local oscillator and signal beams is 



MEMS: microelectromechanical system 
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crucial. For this reason, identical collimators in both beams are employed. 
Careful alignment then allowed to achieve interference contrast up to 90% 



(figure 48). 
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Fig. 48. Interference contrast of the homodyne receiver. "Signal" and "local oscilla- 
tor" traces are recorded by manually rotating the last A/2-waveplate in front of the 
last polarizing beam splitter cube in steps of 10°, with the other input beam blocked. 
For this test, the powers in both beams are adjusted to yield a maximum signal of 
400 mV. If both beams are opened, and the detection unit is properly aligned, a 
total signal amplitude of 1440 mV, corresponding to 90% interference contrast, is 
recorded. 

The phase of the local oscillator is locked by using the DC-part of the homo- 
dyne signal from the balanced receiver, forcing the average detected differential 
photocurrent to zero for a phase measurement. Note that this is slightly dif- 
ferent to the technique described in |4.2.H where an orthogonal polarization 
component (in both signal and LO beams) was used to derive the LO phase 
error signal. In this experiment, the length of the employed fibers renders this 
approach problematic due to apparent drifts in the polarization rotation in 
the fibers. Note also that it is necessary to match the total length of the fibers 
in both arms of the interferometer (excluding the potential propagation length 
in the toroid) in order to measure the correct signal quadrature at all Fourier 
frequencies. 



Figure [49] shows typical displacement spectra of the 65.3 MHz- RBM of a silica 
microtoroid probed by homodyne spectroscopy with the Ti:S laser. The sample 
is held in the cryostat, at a temperature of 2.4 K, corresponding to (n) ^ 770. 
Sensitivity at the 10"-^^ m/v/Hz-level, as well as a SBR of nearly 20 dB (in 
noise power) is achieved, in spite of a mechanical quality factor which was as 
low as 540 in this case. Note that displacement noise measured at resonance 
is the sum of the thermal noise in the mechanical mode plus the imprecision 
background noise due to quantum noise in the detection, which is reduced at 
higher probing powers. 



7.2 Cooling towards the quantum ground state 



To demonstrate the performance of combined optical and cryogenic cooling, 
we used a 52 /xm-diameter sample with a WGM of 5.5 MHz intrinsic decay 
rate and 9 MHz mode splitting, loaded down to /^/27r ^ 19 MHz using the 
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Fig. 49. Displacement noise spectra of the RBM of a silica microtoroid held at a 
temperature of 2.4 K. The power used to probe the RBM was about (red 
trace), 10 fiW (green trace) and 100 /xW (blue trace). Points are measured data, 
lines are Lorentzian fits. Figure from ref. ([Schliesser et al. (2009)). 



fiber taper, corresponding to a finesse of about 70,000. A room-temperature 
reference measurement (figure 49) of its RBM at 65.2 MHz reveals an effective 
mass of 5.6 ng, in good agreement with the simulated value of 4.9 ng. Note 
that the experimentally determined mass can only be as accurate as is the 
displacement calibration and the knowledge of the actual temperature of the 
RBM. We generally estimate the resulting systematic error to a level of 30%, 
arising from imperfections in the modulation scheme and temperature changes 
induced by dynamical backaction or absorption-induced heating. 



Using this effective mass, it is possible to derive the mode temperature from 
subsequent displacement measurements during cooldown of the cryostat. The 



results again confirm proper thermalization of the sample (figure 49), down 
to an occupation of (n) ^ 550. For further cooling, the Ti:S laser is tuned 
to the lower sideband at A = — f^m- The signal from the homodyne receiver 
can still be used to measure mechanical displacements, however, due to the 
detuned operation, the spectrum of the phase quadrature now displays a more 
complicated dependence from the mechanical noise spectrum. In contrast to 



the simple expression (|118|) for A = 0, we now obtain 

sTi^) = 1 



+ 2A2(277e - l){K/2f + (277e - l)2(«:/2)2((«:/2)2 + Q^) 



SxX (^) 5 



A4 + 2A2((«;/2)2 - Q2) + ((^/2)2 + Q2)2 

(180) 

where the first term, equal to unity, again is due to quantum noise. Impor- 
tantly, the transduction of a frequency modulation into Sqq has the same 
dependence on the relevant experimental parameters as the transduction of 
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Fig. 50. Left panel: Mode temperature and occupation number of the RBM, mea- 
sured using noise thermometry, versus the temperature of the buffer gas in the 
cryostat, including a reference measurement at room temperature (red point). Green 
data points are taken with the laser tuned to the optical resonance, so no optical 
cooling takes place. When tuned to the lower mechanical sideband (inset), addi- 
tional laser cooling is observed (blue points). Panels on the right show a selection of 
displacement noise spectra, from which the mode temperatures of the measurements 
were derived. Apart from the calibration peak these spectra also reveal a nearby 
second mechanical mode. Nonetheless, the spectrum of the RBM can be well-fit 



with a Lorentzian (thin lies). Figure from ref. (Schliesser et al. ( 2009| )). 



displacements Sx 
detuned case. 



so that the calibration scheme remains valid also in the 



Rewriting equation (180) as a quantum-noise induced imprecision in the dis- 



placement measurement one obtains 
A2 + (277e-l)2(«:/2)2 



Aa?gQr}cK 

A4 + 2A2((«:/2)2 



^2) + ((/t/2)2 + r22)2 



A4 + 2A2(2r],-l)(«:/2)2 + (277e 



l)2(K/2)2((«:/2)2 + r22)- 

(181) 



For K ^ Feff , which is typically in the case in the experiments presented here, 
the spectral shape of this function can be assumed to be flat over the frequency 
range of interest, and we find as a useful approximation 



Qim,qn/Q 



(182) 



in the resolved sideband regime (A = — f^m S> /^). It is noteworthy that this 
value is only a factor of 4 higher than expected for resonant probing A = 
in this regime. However, in the resolved sideband regime, only a fraction of 
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{k/AD.^Y' of the launched power is coupled to the cavity. 



Figure [50] shows displacement spectra and corresponding mode temperatures 
recorded during such a cooling run. As expected, the damping of the mechan- 
ical modes is optically increased, when the laser is detuned to A = — f^m- 
Note that for these experiments, active stabilization of the laser frequency to 
the optical sideband of the WGM has not been necessary. Instead, a second 
electronic spectrum analyzer was used to demodulate the homodyne signal at 
the frequency f^mod of the calibration modulation. As expected from equation 



(181), this signal displays local maxima at A = O,±f^mod- After the laser is 
tuned to A = — f^m, the system is sufficiently stable during the averaging of 
typically several tens of seconds (small drifts can be manually compensated). 

For the highest launched powers of ^ 0.2 mW, the total damping rate was 
increased to reff/27r = 370 kHz. At the same time, the mode temperature was 
reduced to = 200 ± 60 mK, corresponding to an occupation number of 
(n) = 63 ± 20 dSchhesser et~aL] ( [20091 )). The error interval is due to potential 



calibration errors in the displacement measurement and the determination of 
the effective mass. This value is comparable to the lowest occupation number 
achieved in nanoelectromechanical systems, yet those experiments have to rely 
more heavily on advanced cryogenic machinery (Naik et al. ( 2006[ )), and suffer 
from the insufficient signal-to-noise ratio. 



7.3 Assessing measurement backaction 



The strong suppression of thermal noise in the measurements presented here 
allows an interesting analysis of the data from the perspective of quantum 
measurement (Braginsky and Khalili ( 1992D ). As discussed in section 4.1, fun- 
damentally, the achieved total uncertainty in the measurement of the dis- 
placement of the mechanical oscillator is limited by two effects: imprecision 
and backaction. The imprecision in our experiment is given by detection shot 



noise, in the detuned case discussed here it is described by equation (181). 



On the other hand, the light used to measure the mechanical oscillator exerts 
backaction on it. Inevitably, intracavity photon number fluctuations give rise 



to a fluctuating force, the spectrum of which can be derived as (Schliesser 
eFaL] ( [2009( )) 



(183) 



^ A4 + 2A2((/^/2)2 - n^) + ((/^/2)2 + n^y ' 

These force fluctuations are usually referred to as quantum backaction, sim- 
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plifying to 



Qba,qn 



2^0-Pin??c^ 



(184) 



in the resolved-sideband regime (|A| = ^ /^). Note that the spectra of 
imprecision and backaction noise, as weU as their possible correlation Sxf{^) 
reflect properties of the measurement device (the cavity pumped by a laser 
field) — independent of the mechanical oscillator. 



Other sources of measurement backaction include excess noise in the intra- 
cavity photon number, which may arise from laser frequency fluctuations 



dSchhesser etaL] ( |2008b[ ); [Di6si| ( |2008| ); |Rabl et al.| ( |2009| )). This effect can 
be ruled out in this work as the employed Ti:S is known to exhibit only quan- 
tum fluctuations at the Fourier frequencies of interest. Another possible source 
of measurement backaction is heating of the torus due to light absorption. This 
increases the temperature of the structure, and raises the level of Langevin 
force fluctuations driving the mechanical oscillator. 



A series of cooling measurements (flgure 51) indeed reveals a deviation from 



the simple relation (147), which can however be reproduced by introducing a 
heating term 

r„(r) 



with = T + ATabs, where ATabs is proportional to the power circulating 
in the cavity. In addition, we have taken also the temperature dependence of 
the mechanical damping r^{T^) into account. For the 65 MHz-oscillator below 
2K, with ( |Arcizet et al] ( p09| ) dV^/dT ^ 27rl6kHz/K, we flnd a heating of 
about 5 K per Watt of circulating power. Similar values were extracted from 
studies of optical multi-stability at low temperature at a wavelength of 1.5 /xm, 
corroborating the attribution of the observed backaction effect to laser-induced 
heating. 

A comparison with a second cooling run with a different sample with Vt^/2T{ = 
121.7 MHz and = 2,200, but a broader WGM resonance (/^/27r = 155 MHz) 
emphasizes the importance of the resolved-sideband regime for the efficiency 
of cooling in the presence of laser absorption: A signiflcantly more pronounced 
heating effect prevents reaching occupation numbers below (n) = 100, in spite 



of the higher mechanical frequency (flgure 51) 



From the data of the (n) ^ 63-cooling run we can now extract quanti- 
tative values of the corresponding backaction force fluctuations. As a very 
conservative upper limit, we may use the total thermal force fluctuations 
Sf§^{n) = 2meff^^A;B^^ and flnd a value of ^Sf§^{n^) = SfN/^Ife. If we 
consider only the temperature rise ATabs caused by the laser absorption as 
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Fig. 51. Resolved-sideband cooling of two samples with frequencies of 
f^m/27r = 65.1MHz (left panel) and f^m/27r = 121.7MHz (right panel). The graphs 
show the mechanical linewidth Feff (abscissa) versus the derived mode tempera- 
ture Tm (ordinate) during a cooling run. Open points correspond to measurements 
with the laser tuned close to the optical resonance (no optical cooling), and filled 
points to measurements with the laser tuned close to the lower sideband. When 
varying the power of the cooling laser, both linewidth and mode temperature are 
changed. Deviation from the simple cooling behavior (blue line) is attributed to an 
increase of the structures' temperature (red dashed line), taken into account in a 
more elaborate model (green dashed line). Scatter is due to varying operation con- 
ditions, uncertainty in phonon occupation for each point is < 30%. Figure from ref. 
dSchliesser et al.| ( |2009D ). 



the backaction of the measurement, a lower value of y Sp%^{Qrn) = 4 fN/v^Hz is 
found. In these assessments, we benefit from the low occupation which allows 
us to extract the effect of measurement backaction, as it is large enough to be 
observed on top of the background of the thermal noise. 



It is interesting to compare these findings with fundamental limits. Quantum 
mechanics imposes an inequality on imprecision and backaction noise, which, 
for the particular case of an optical measurement of a mechanical oscillator's 



displacement, can be written as (Braginsky and Khalili ( |1992| )) 



r/im 



4 ' 



(186) 



This relation can be considered a manifestation of the Heisenberg uncertainty 



principle in the context of continuous position measurement (Braginsky and 
Khalilil ( [19921 )). 



Taking the force noise extracted from our data, and the experimental impre- 

1.4am/v^Hz achieved in the same measurement, an up- 



cision of ^JS'^{Qrn) 
per limit from the backaction-imprecision product of 
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220/1/2 is found. Considering only the absorption-induced heating as a back- 
action mechanism, an even lower value of 



(n^) 'S^%^{n^) ^ 100 /i/2 is 

found from our experiments. This is an order of magnitude lower than the val- 
ues achieved with nanomechanical oscillators cooled in dilution refrigerators: 
Readout with an atomic point contact (Flowers- Jacobs et al. ( 2007] )) achieved 
a backaction-imprecision product of 1700±400 /i/2, while measurements using 
a superconducting single-electron transistor (Naik et al. ( 2006D ) have achieved 
a value 1^ of - 800^/2. 



Note that the number quoted in this manuscript, 15 /i/2, is a theoretical estimate 
if shot-noise limited detection was possible. As discussed in the supplementary in- 



formation, the actual imprecision noise in the experiment was 50-times higher (Naik 



et al. 



(2006); Clerk et al. (2008)). 
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8 Conclusion 



In this chapter, we have reported on optomechanical interactions in high- 
finesse optical whispering gallery mode resonators. Careful analysis and un- 
derstanding of these devices' properties has allowed us to optimize their perfor- 
mance for the purpose of cavity optomechanics. Among the various systems 
now designed and studied in this context, they offer a unique combination 
of high-frequency (30-120 MHz), high-quality (Qm up to 80,000) mechanical 
modes and ultra-high finesse (intrinsic finesse J^o ^ 0.9 • 10^) optical reso- 
nances. These key figures rival even the best optical cavities developed in the 
context of cavity quantum electrodynamics (Rempe et al. (19921)) in terms of 
finesse, and the mechanical quality factors of state-of-the art nano- and mi- 
cromechanical oscillators (Naik et al. (2006); Verbridge et al. ( 2007] )) in the 
same frequency range. 



Together with an intrinsically strong parametric coupling of optical and me- 
chanical degrees of freedom — with optical resonance frequency shifts of typ- 
ically more than 10 kHz per femtometer displacement — this has enabled ul- 
trasensitive optical monitoring of the mechanical mode. Adapting powerful 
tools from quantum optics as optical displacement meters limited only by op- 
tical quantum noise, sensitivities at the level of 10~^^m/v/Hz were achieved 
(Schliesser et al. ( 2008b|a| )). This experimentally demonstrated imprecision 
is below the expected noise level associated with the zero-point fiuctuations 
of the mechanical mode, which so far has been achieved only with much 



lighter nanomechanical oscillators ( |Naik et al. (2006); Clerk et al. ( 2008[ )) 



Reaching such a sensitivity is a crucial precondition for the experimental 
verification of the concepts of quantum measurement in the context of dis- 
placement measurements, such as quantum backaction and the emergence of 



the standard quantum limit (Braginsky and Khalili (1992); Tittonen et al. 



(1999); Schwab and Roukes (2005)). Furthermore, this exquisite sensitivity 
(Schliesser et al. (2008a)) may also be exploited to monitor nanomechanical 
oscillators brought into the near-field of the cavity mode, such as silicon ni- 



tride nanobeams or -membranes (Anetsberger et al. (2009)), graphene sheets 
or diamond nanowires. Beyond mechanical effects, this sensitivity has allowed 
us to study fundamental thermal noise mechanisms, such as thermorefractive 
noise, which are of interest for the application of silica microresonators as fre- 



quency references (Vassiliev et al. (1998); Matsko et al. (2007)), as biophysical 
sensors ( [Armani et al.| ( |2007D ; [Schroter et al.| ( |2008D ; |Vollmer et al.| (|2009^) or 
for the proposed demonstration of Kerr squeezing (Kippenberg et al. ( 2004[ )). 



Due to the high-finesse, the dramatically enhanced intracavity radiation pres- 
sure exerts a readily detectable force on the mechanical mode. This effect was 
directly measured in a "pump-probe" -type measurement, in which the dis- 
placements induced by modulated radiation-pressure is probed with a second 
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laser. It was shown that the nonhnear cross-couphng of the two hght fields is 
strongly dominated by radiation pressure induced mechanical displacement, 
and more than two orders of magnitude stronger than the well-known Kerr 
effect (Schliesser et al. ( 2006D ). Furthermore, radiation pressure has also been 
shown to induce a modification of the dynamics of the mechanical mode, 
changing both its effective spring constant (optical spring) and its damping 



(amplification and cooling). Predicted as early as 1967 by Braginsky (Bra- 



ginskii and Manukin (1967)), this dynamical backaction was systematically 



studied over a wide regime of experimental parameters (detuning, photon stor- 
age time, mechanical oscillation period). Our experiments have demonstrated, 
for the first time, efficient optical cooling of a mechanical mode induced by 
dynamical backaction, both in the "Doppler" (Schliesser et al. (2006)) and 
the resolved-sideband regime (Schliesser et al. ( 2008b[ )). These techniques are 
now widely employed in experiments which aim to demonstrate ground-state 
cooling of a mechanical device (figure 52). 



In our experiments, we have identified several important barriers on the way 
towards this ultimate goal. Frequency noise of the driving electromagnetic 



field, practically relevant in many systems (Schliesser et al. (2008b); Diosi 



(2008); Rabl et al. ( 2009] )), is shown here to be eliminated by resorting to a 



quantum-noise limited laser system. Furthermore, much in contrast to trapped 
atoms or ions, even very high-Q macroscopic oscillators are not very well 
isolated from their finite-temperature environment. Fluctuating thermal forces 
thus compete with laser cooling, limiting the occupancies achieved in room- 
temperature experiments to a few thousand quanta (figure 52). 



To overcome this limit, we have implemented resolved-sideband laser cooling in 
a cryogenic environment (Schliesser et al. (2009)). The mechanical oscillator is 
cooled to an occupation number of (n) ^ 63 ±20. The optical detection scheme 
provides small enough displacement imprecision to enable monitoring even 
such an ultracold oscillator with appreciable signal-to-noise ratio. Additionally, 
due to the low occupation number, we are able to extract an upper limit on 
the backaction of the displacement measurement. The product of backaction 
and imprecision noise lies only a factor of 100 above the fundamental quantum 



limit (Braginsky and Khalili (1992)), and constitutes the lowest value reported 



in the literature (Clerk et al. 



p08| )). 



In conclusion, we have introduced an optomechanical system suited for the 
exploration of quantum effects in mesoscopic mechanical oscillators, by ex- 
perimentally demonstrating (i) the possibility to monitor its displacements 
with an imprecision below the standard quantum limit, (ii) di combination of 
cooling techniques based on radiation-pressure and standard cryogenics ca- 
pable of cooling the oscillator close to its quantum ground state, and (Hi) a 
near-ideal operation of the displacement transducer in the sense of quantum 
measurement (imprecision-backaction product). Leveraging quantum optical 
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Fig. 52. Cooling experiments performed in different laboratories around the globe, 
including both experiments based on dynamical backaction and active feedback 
cooling. The cooled oscillators span about 8 orders of magnitude in frequency. Open 
symbols indicate the reservoir temperature of the experiments, distinguishing cryo- 
genic from room-temperature experiments. Full symbols indicate the lowest mode 
temperature achieved when optomechnical cooling is applied. The dashed line indi- 
cates the temperature for which hQ^n ^ ksT. The individual results are described 
in references: AURI GA ([Vinante et aT] (|2008D), C ornell ( [Hertzberg et al.| (2009^), 



MIT dCorbitt et aL| ( |2007a|bD ), MPQ (,Schliesser et al. 



dBrown et al.| ( |2007D ), O regon ([Park and Wang| (|2009|), Sta nford ( jWeld and Ka- 



IBM (IPoggio et al.| (|2007D), IQOQI ( |Gigan et al.| ([2006D; problacher et al.| (|2008l 
2Q09b )), JILA dTeufel etaTI ([20081)), LKB ([Cohadon et al.l (|1999D; [Arcizet et al. 
( |2006bpi|)), LMU dHohberger Metzger a nd Karrai| ([2004^7 F'avero et al.| ( |2007D )7 



(2006, 2008b 



2009)), NIST 



pitulnik] d20"Q6l ) ), UCSB ([Kleckner and Bouwmeester] d2006D), U MD dNaik et al. 
d2006D ), UWA dMow-Lowry et al.| d2008D ), Yale dThomson et al.| ( |2008l )). 
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techniques, we have therefore established a route into the quantum realm of 
mesoscopic oscillators, which previously was thought to be uniquely accessible 



with nano-electromechanical systems (Schwab and Roukes (2005)). 
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9 Outlook 



Notwithstanding such progress, true quantum effects of radiation pressure and 
the quantum states of a mechanical osciUator have not been observed in any 
experiment today. In our experiments involving WGM micoresonators, it 
appears that cooling to the quantum ground state, and the observation of 
quantum backaction are closest to fruition. In the following, we briefly outline 
the strategies pursued in our laboratory to achieve these goals, and, concluding 
this chapter, point to other interesting optomechanical eflFects. 

The most severe antagonist to ground state cooling is the intrinsic damping 
of the mechanical mode, coupling it to the environment at a finite temperature 
T. This coupling tries to maintain the mode and the environment in thermal 
equilibrium, essentially feeding the power T^kB{T — T^) into the mechanical 
mode, which has to be continuously removed by laser cooling. There are two 
obvious ways to reduce this load: reducing the environment temperature T, 
and reducing the mechanical damping T^. While the former is a technical task, 
the latter exhibits complex dependence on material, geometry and operating 
conditions. 



In particular we have observed (sections 3.2.4 and 7.2) a strong increase of the 



mechanical damping at cryogenic temperatures due to relaxation of two-level 



systems (TLS) present in silica material (|Pohl et al. (2002); Arcizet et al. 



(2009j)). This fiaw is known, however, to ameliorate at suflBciently low temper- 



atures < 1 K, at which the damping due to TLS falls oflF with oc Q^/T^. 
Thus, lowering the operating temperature of the experiment may be expected 
to enhance the cooling performance with a scaling up to T~^. For this reason, a 
^He-cryogenic system is presently being tested in our laboratory. A base tem- 
perature of 600mKand significantly higher mechanical quality factors (> 10^) 
than in the ^He-cryostat were already measured. Resolved-sideband cooling, 
with the cooling rates of 1.5 MHz already demonstrated, would place a 70 MHz 
RBM directly in the quantum ground state. A very crucial issue in this con- 
text, however, is the suppression of heating by laser absorption, requiring very 
high-quality samples and operation deeply in the resolved sideband regime. 

Optical absorption could be reduced by using ultrapure crystalline materi- 
als for the WGM resonators. Indeed, pioneering work at the Jet Propulsion 
Laboratory has resulted in WGM resonators made of quartz (?), MgF2 and 



CaF2 (Grudinin et al. (2006)) having optical quality factors up to 6- 10^^. Ma- 
chining and polishing WGM resonators down to a diameter of 80 fim has also 
been achieved. At the same time, the pristine crystalline structure avoids me- 



There are, however, experiments that could be considered "quantum simulations" 



of mesoscopic mechanical oscillators — based on ultracold atoms ( |Murch et al.| ( |2008| ); 



Brennecke et al. 



(2008)) 
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Fig. 53. Whispering gallery mode resonator made of calcium fluoride. The rims of 
the protrusions support WGM with extreme photon storage times. Left panel is 
a photograph of a 6mm-diameter prototype resonator produced at MPQ. Right 
panel shows a ringdown trace for a critically coupled resonator. The ringdown time 



of 2.5 /xs corresponds to an intrinsic quality factor of 1.2 • 10^ (Hofer et al. (2009)). 



chanical losses due to two-level systems. At ^ 100 kHz- frequencies, mechanical 
quality factors Qm ^ 10^ were measured in bulk CaF2 samples, both at room 



and cryogenic temperatures (Nawrodt et al. (2007)) 



We have therefore started to explore the optomechanical properties of crys- 
talline WGM resonators in our laboratory. Several diflFerent geometries were 
fabricated using a precision lathe and polished using diamond slurry. An exam- 
ple of such a resonator is shown in figure [53) We have achieved optical quality 
factors up to Qq = 1.2-10^°, corresponding to a linewidth of /^/27r = 24 kHz in a 
R = 1.8mm resonator (intrinsic Finesse J^q = 760,000). We have also produced 
an 800 /xm-diameter, 100 /xm thick disk and achieved Qq ^ 10^ (J^o = 400,000). 
Optical transduction techniques described in section [4] are used to measure 
mechanical modes in these structures, which are typically found in the range 
between 0.5 and 5 MHz. The highest measured mechanical quality factor was 
Qm = 136,000. Already an encouraging value, there is strong evidence for this 
value to be still limited by clamping losses which can be mitigated by a more 
suitable design and suspension of the resonators. Furthermore, the rather high 
effective masses (around 600 /ig) and the slightly weaker optomechanical cou- 
pling {\go\/27T ^ 1.5kHz/fm) call for further miniaturization of the structures. 
For example, a 80 /xm-diameter, 10 /xm thick disk would possess m^ff = 90 ng 
and Q^/27T = 63 MHz. If Qm ^ 10^ can be reached, the power required to cool 
such a device from T = = 1.6K to ^ HQ^/kB is as low as 10 /xW in 
the resolved-sideband regime. Heating due to absorbtion is likely to be totally 
negligible considering the optical quality of the crystals. 

We finally point out that many other crystalline materials may be amenable 
to this approach, combining optomechanical coupling with yet other function- 



105 



ality. As an example, we have tested various polished diamond spheres, and 
have observed WGMs with quality factors up to 3 • 10^ in a 3 mm-diameter 
sphere. 

A somewhat opposite approach consists in decoupling optical and mechanical 
degrees of freedom. Keeping the silica WGM resonators as optical cavities, 
but placing an external mechanical oscillator in the near-field of the WGM 
allows for independently engineered mechanical properties. For example, light- 
weight (meff ^ 1 pg) nanomechanical oscillators such as strained SiN strings, 
have been shown to result in Qm ^ 10^ at ^ 1 MHz-resonance frequencies 



(Verbridge et al. ( 2008| )). Placed in the near-field of a WGM, movement of 
the polarizable oscillator in the optical field gradient induces frequency shifts 
on the order of l^'ol ^ lOMHz/nm, without inducing detectable optical loss 
of the WGM. A measurement imprecision below the noise level associated 
with zero-point fiuctuations and dynamical backaction induced by the optical 
dipole force have been demonstrated with such a system in our laboratory 
dAnetsberger et al.| ( |2009| )). 



These results are very encouraging for studies pertaining to the observa- 
tion of quantum backaction. Usually masked by the much stronger thermal 
noise, this effect may become detectable if the ratio 5'^^^^(f^m)/5')^^F(^m) ^ 
a^g^nli?^'^ /2mQ^V^kT (assuming the limit ^ ^/2) approaches unity. 
Leveraging the low mass and high quality factor available with nanomechan- 
ical oscillators, it appears feasible to approach this regime, with moderate 
optical probing powers of ^ 100 /xW even at room temperature. To differenti- 
ate the added noise induced by quantum backaction from potentially present 
absorption induced heating, correlation measurements ( Verlot et al. ( 2008D ) or 
backaction cancellation (Caniard et al. ( 2007b| )) between two oscillators with 
slightly different frequencies may be employed. 

A final example of near-future research projects enabled by the progress pre- 
sented in this thesis is an investigation of the regime of strong optomechanical 
coupling ("Dobrindt et al. ( 2008] )) in the yet unexplored optical domain. As our 
analysis shows, the strong — but tunable — coupling of the optical mode to the 
mechanical mode via the mean photon field a gives rise to features very simi- 
lar to the effect of electromagnetically induced transparency (EIT) in atomic 

( [205^ ). 



physics (Zhang et al. (2003); Fleischhauer et al 



The basic idea of such an experiment is illustrated in figure [54} The strong 
field a of the coupling (formerly "cooling" ) laser oscillates at frequency (jJ\ = 
cUc — ^im and thereby couples the levels 2 ^ 3 by processes in which a phonon is 
removed upon the addition of a photon to the intracavity field (red-sideband 
transitions). A second, very weak laser oscillating at uj^ = uj\ + f], probes 
"carrier" transitions 1 2, in which originally no phonons are added or 
removed. Due to the strong optomechanical coupling induced by the coupling 
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Optical frequency 



Fig. 54. Optomechanically induced transparency. Left: A coupling laser of frequency 
cj\ drives red sideband transitions between (sets of) levels 2 and 3, which differ both 
in the number of optical and mechanical excitation quanta (each time a photon is 
added to the coupling field, a phonon in the mechanical degree of freedom is re- 
moved). A probing laser tuned close to the 1^3 transition of the unperturbed sys- 
tems actually interacts with a superposition of levels 2 and 3, leading to the opening 
of a transmission window in the case of an optomechanical system. Right panel shows 
the resulting power transmission of the probe laser through a taper-toroid system 
for A = —^m versus the detuning {^ — ^^)/27r of the probe laser. Numerical param- 
eters are typical for silica microtoroids, with fim/27r = 40 MHz, rm/27r = 1.3 kHz, 
i^/27v 5 MHz, R =, meff 10 ng, i? = 40/xm and s?^h(jj\ 300 /xW (red line) or 
s?^hijj\ = 0/iW (green line). 



field, it effectively interacts with a hybrid optomechanical resonance. Similar 
to EIT, this hybridization opens up a tunable transmission window for the 
probing laser at the center of the optical resonance. We therefore refer to this 
effect as "optomechanically induced transparency" . 



A simple model for this scheme can be set up directly from the Langevin 



equations (50)-(52). Neglecting thermal and quantum noise for an elementary 



analysis, one obtains for the (power) transmission of the probing laser 



1 - 



-i{A + Q) + i^/2 + 2Af{Q) 



(187) 



with 



i{/\-n) + K/2' 



(188) 



At a basic level, this effect can be understood from the fact that for cjp — 
uj\ = Q ^ r^m, the beat of coupling and probing laser drives the mechanical 
oscillator. The anti-Stokes field generated, in turn, from the coupling field 
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interferes with the incoming probing hght. 



The formal analogy with atomic EIT is even more obvious in the limiting case 
of: (i) Si high-(5 oscillator ^ (H) the resolved-sideband regime ^ k 
and (Hi) the detuning A = — f^m- Using the abbreviation A' = — one 
can then simplify the equation system to 



hgod 



(189) 
(190) 



directly analogous to the EIT case (Milonni (2005)), with the well-known 
solution 

^^(^) = '^^'^ ....... • (191) 



(- 



Mechanical and optical oscillator play the role of two dipole transitions, where 
a pump laser couples two of the involved three levels. As in EIT, the modifi- 
cation of the transmission for the probe laser may be used for tuning optical 



group velocities to generate slow or fast light (Kash et al. (1999); Hau et al. 



(1999); Milonni (2005)), or, in designs hosting large arrays of optomechanical 



systems, even storage of light in mechanical states (iFleischhauer and Lukin 



(2000); Liu et al. (2001); PhiUips et al. (2001)) may be envisioned. 
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